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' Abstract 

In this paper we describe the structure of a class of two-component scalar field models in 
a (1+1) Minkowskian space-time which generalize the well-known Montonen-Sarker-Trullinger- 
Bishop -hence MSTB- model. This class includes all the field models whose static field equations 
are equivalent to the Newton equations of two-dimensional type I Liouville mechanical systems 
qq [ with a discrete set of instability points. We offer a systematic procedure to characterize these 

£S| ■ models and to identify the solitary wave or kink solutions as homoclinic or heteroclinic trajecto- 

ries in the analogous mechanical system. This procedure is applied to a one-parametric family 
■ of generalized MSTB models with a degree-eight polynomial as potential energy density. 

1 Introduction 

Over the last decades topological defects or solitary waves have played an essential role in the 
\ explanation of new phenomena in diverse branches of Science, e.g., Cosmology and Condensed 
Matter Physics. These non-linear waves are non-dispersive spatially localized solutions of nondinear 
field equations describing a given physical system. Focusing on one-dimensional wave phenomena, 
we shall deal with variations of the relativistic non-linear Klein-Gordon PDE [lj 
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for N scalar fields (p a (x,t). The non-linear nature of these equations is the key point to finding 
unexpected a priori non-dispersive wave solutions. The simplest examples of relativistic topological 
defects of this kind are the soliton of sine-Gordon theory and the kink of the ((^)| model. These 
systems encompass only one scalar field and provide theoretical support to explain, for instance, 
the appearance of superconductivity in type II materials [21 [3], electric charge fractionization in 
trans-polyacetylene (CH) X [3], the Josephson effect, [5], etcetera. In several branches of Physics 
ranging from condensed matter physics to cosmology, from fluid dynamics to relativistic quantum 
field theory, more complex systems involving several scalar fields arise. The search for topological 
defects, solitary waves or kinks in models which involve several scalar fields, as Rajaraman notices 
[6], is a very difficult endeavor: This already brings us to the stage where no general methods are 
available for obtaining all localized static solutions (kinks), given the field equations. However, 
some solutions, but by no means all, can be obtained for a class of such Lagrangians using a little 
trial and error. The move from one-component scalar field models to two-component scalar field 
models is thus a very important qualitative step, passing from the integration of a single partial 
differential equation to the integration of a system of partial differential equations. Although there 
are no general methods to tackle the problem of solving systems of non-linear partial differential 



1 



equations, some strategies for finding soliton solutions have been developed in both the physical 
and mathematical literature over the last thirty years. 

We now point out two of these strategies that lie at the core of the conceptual framework to be 
developed in this paper: 

• Since the models we are interested in exhibit Lorentz invariance, see (pQ), the search for kinks 
or traveling waves is tantamount to the solution of an analogous mechanical problem where 
the unit-mass point particle motion in a certain potential V(fa, fa) must be identified: this 
analogy is established by thinking of the variable x as "time" whereas the field components 
(01) fa) are reinterpreted as the coordinates of the particle [Bj moving in the Euclidean plane 
under the influence of the potential V{fa,fa) = —U(fa, fa). Moreover, the field theoretical 
energy becomes the mechanical action. This identification means that finite energy field 
configurations correspond to finite action (in infinite time) paths in the analogous mechanical 
system. 

By this token one must solve only a system of ordinary differential equations, instead of 
partial differential equations. Traveling wave solutions to the PDE system are thus obtained 
by applying Lorentz transformations (or Galilean, in non-relativistic systems) to the static 
solutions of the ODE system. 

• In some models the potential density energy can be written as half the square of the norm of 
the gradient of a "superpotential" W(fa,fa): 

TJI1 1 , 1 (dW dW dW dW 

U ^ ) = 2[dfa'dfa + dfa'dfa 

Knowledge of such a function, W(fa,fa), allows the Bogomolny'i-Prasad-Sommerfield ar- 
rangement [7] for the energy functional: 

pu 1 f, s^\dfa dW] 2 f ^dWdfa 

Because it is an exact differential if W is well behaved, the second term 
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where by K.(x) £ Maps(K,M. 2 ) we denote a non-singular curve in M 2 , in ^ only depends on 
the values of W at the end points of the curve: W[fa(±oc,t),fa(±oo,t)]. If finite energy is 
required, the following conditions at spatial infinity, (M, must be satisfied: 

{fa(±oo,t),fa(±oo,t))eM , ^(±oo,t) = = ^(±oo,t) 

where Ai is the set of zeroes of U, assuming a non-negative U. Because the temporal evolution 
is a homotopy transformation, when A4 is a discrete set W[0i(±oo, t), fa(ztoo, t)} is time- 
independent, and T is a "topological charge". 

The second term in ([3|) is semi-definite positive and the absolute minima of the energy are 
solutions of the first-order system of ordinary differential equations: 

dfa = 9W a = l2 

dx dfa 
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Non-dispersive extended solutions of (|1J) are also stable solutions of the second-order ODE 
system §5§ obeying the static field equations: 

<Ha = dU_ 

dx 2 dcj) a ' 1 j 

usually referred to as solitary waves, kinks, or one-dimensional topological defects. 

The link between these two strategies is hidden in equation This equation is no more 
than the reduced Hamilton- Jacobi equation for zero mechanical energy trajectories of the 
analogous mechanical system. As a consequence W is the Hamilton characteristic function, 
whereas the term superpotential comes from the fact that only models with potential energy 
of the form written in ([2]) are susceptible to super symmetric extensions. 

In this way, the search for solitary waves in N scalar field theory is tantamount to the 
solving of a iV-dimensional mechanical system. Thus, the most favorable situation is to 
deal with integrable mechanical systems. If N = 2, there is a lot of information about 
classes of integrable mechanical systems, see e.g. [41] . In particular, we shall choose field 
theoretical models such that their analogous mechanical system is of Liouville Type I: those 
two-dimensional mechanical systems such that their Hamilton-Jacobi equation is separable 
using elliptic coordinates. The simplest and best studied (1+1) dimensional two-component 
scalar field theory model of this type is the MSTB model (after Montonen, Sarker, Trullinger 
and Bishop, who first addressed this theory in \17\ I18j). 

Our goal in this paper is to analyze the structure of this broad class of systems of the MSTB 
(the tip of the iceberg) type as well as to study their manifold of solitary waves. Applications 
of multi-component kinks to describe interesting physical systems abound in the literature, 
see [7-15]. 



1.1 The MSTB Model 

The MSTB model is a physical system with a proud history. In 1976 Montonen [IT], searching for 
charged solitons in a model with one complex and one real scalar field, discovered, by fixing the 
time-dependent phase for the complex field, a (l+l)-dimensional two real scalar field theory that 
provided the basic neutral solitons. When such basic neutral solitons are embedded in the bigger 
system charged solitons arises. Thus, the MSTB model is a (l+l)-dimensional relativistic theory 
of two scalar fields with potential energy densit}Q: 

i n 2 

U MS TB[Mx,t),Mx,t)] = -((P 2 (x,t) + ^ 2 (x,t)-l) 2 + —<j) 2 (x,t) . (6) 

Considered as a function of the two scalar fields 0i and (p2, C^mstb is a fourth-order polynomial 
isotropic in quartic but anisotropic in quadratic terms. There is only one non-dimensional param- 
eter 0, which determines the intensity of the anisotropy (as well as the time dependence of the 
phase of the complex scalar field) . 

The kink solutions asymptotically connect elements of the set of zeroes of £7mstb(^i> 4>2), which 
in this case are two : M. = {(4>f = 1, c/>^~ = 0); {4>i = — lj 4>2 = 0)1- ^ n ^ s seminal paper [17] , the 
author identified two solitary waves or kinks, joining the points <^>+ and <j>- asymptotically in the 
range of the coupling constant: £l 2 G (0, 1). The first kink was the old kink of the §\ theory with 
only one real scalar field: 

^ K1 (x)=tanh^ , ^ K1 W = ; x= ~ ^ ~ (7) 

Vl — v z 

1 We shall use non-dimensional coordinates, fields, and parameters throughout the paper. 
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where a = 0, 1, v is a velocity slower than the speed of light c = 1, and 72 £ M sets the kink 
center, see Figure 1. The dependence on x is due to the breaking of the translational, Lorentz, and 
reflection symmetries of the model by the kink solution. Because the second field is zero for this 
solution, it is referred to with the abbreviated name of TK1 kink (one-component topological kink): 
the topological charge is non-zero and only one field component is non-null, whereas the solution 
lives in the abscissa axis of the R 2 field space, tracing the orbit: ((fi K1 £ (—1,1), (f^ 1 = 0), which 
starts and ends at different points of M . 

The second and novel kink connects the points 4> + and 0~ through the semi-elliptic orbits 

0? + ff = l , 2 = l-tt 2 (8) 

m the M 2 field space. Solitary wave solutions of this type read: 

cfi K2 (x) = tanhfiz , <£% K2 {x) = (-1)^0 sech tlx ; = 0,1, (9) 

see Figure 1. The abbreviated name for solitary waves of this type is TK2 kink (two-component 
topological kinks) because the topological charge is not zero and the two field components are 
non-null. It is clear that ([8]) becomes a hyperbola when Q 2 > 1 and TK2 kinks do not exist in this 
range. TK1 kinks, however, are still traveling wave solutions for fi > 1. 

To be precise, in a previous paper [19] Rajaraman and Weinberg had identified the first class 
of these solutions and had described the qualitative behavior of the second type of kinks in a more 
general model. Sarker, Trullinger, and Bishop entered the game by establishing the stability of 
the two kinds of solution from a energetic point of view [18]: if Q 2 £ (0, 1) the TK2 kinks, being 
less energetic, are the stable traveling wave solutions and the TK1 kinks should decay to them. 
If Q 2 > 1, obviously the TK1 kinks are stable. In 1979, further analysis of kink stability in this 
model were performed in References [2D] and [2T]. Study of the spectrum of the second-order small 
fluctuation operator around TK1 and TK2 kinks confirmed that the latter are stable and the former 
unstable. The lowest eigenvalue is positive (respectively, negative) for TK2 (respectively,TKl) kinks 
if Q 2 < 1, and positive for TK1 kinks if SI 2 > 1. 

In the same year, Rajaraman [22] discovered a new kind of kink solution for the particular value 
of the coupling constant f2 = \ whose orbits are the circle trajectories: 

These NTK2 kink solutions read 

0NTK2 ( ) = ± 3 h 2 X 1 ^ TK2 (^) = ** sech £ tanh £ 

riw 2 2 2 r2W 2 2 2 

and asymptotically reach <p + , or, <j)~ both at x = ±00, see Figure 1. Hence, these solitary waves 
are referred to as two-component non-topological kinks, or, NTK2 kinks. 

The discovery of this new type of solitary wave prompted several numerical investigations by 
Subbaswamy and Trullinger; their research focused on numerical explorations of the existence, 
energy, and stability of non-topological kink solutions in a wide range of O [23\ [24] . These authors 
numerically found that there exists a continuous family of non-topological kinks of NTK2 type that 
describe closed orbits, all of them enclosed by the ellipse ([8|) in the field space M 2 if ft 2 £ (0, 1). 
Moreover, they numerically checked an astonishing but simple relationship between kink energies: 
the energy of any of the non-topological kinks is the sum of the energies of the two topological 
kinks: 

£[NTK2] = E[TK1] + £[TK2] . (10) 



4 




Figure 1: 0f profile (red solid lines), 4>2 profile (blue solid lines), and energy density (green dashed lines) 
of the TK1, TK2 and NTK2 kinks in the MSTB model. 



This identity, termed the kink mass sum rule, was regarded as a mysterious property of the 
model (consider the kink energy densities in Figure 1). In 1984, Magyari and Thomas |25| noticed 
a very important feature of the model: by finding a second constant of motion, in involution with 
the mechanical energy, these authors showed that the two-dimensional analogous mechanical system 
is completely integrable. In fact, the analogous mechanical model is a particular case of the Gamier 
integrable system, see |26j . Having two isolated unstable points -the minima of U are the maxima 
of V- of maximum potential energy, the heteroclinic trajectories of finite mechanical action joining 
those points are precisely the topological kinks. Homoclinic trajectories, however, starting from 
and ending at the same maximum correspond to the non-top ological kinks. The intelligence of this 
fact plus the existence of two invariants in involution provided an strong hint about the origin of 
the kink mass sum rule as being due to continuous degeneracy in the mechanical action of all the 
separatrix trajectories sharing the critical values of the two constants of motion. 

In 1985, Ito [27] showed that the analogous mechanical system is not only completely integrable 
but Hamilton-Jacobi separable by using elliptic coordinates. The immediate application by Ito 
of the Hamilton-Jacobi principle allowed him to describe the whole kink orbit and kink profile 
varieties analytically by means of complicated expressions between elliptic variables. In particular, 
he verified the sum rule by analytically checking that all the closed trajectories have the same 
mechanical action. 

Although to the best of our knowledge nobody has been able to invert Ito's implicit expressions, 
we recently succeed in doing this by performing a shrewd reshuffling of the kink orbit equations, 
to find: 

1NTK2 / \ (0-l)(l + e 2 ( 1+n )( £ + n T' 1 )) + (fi + l)(e 2n ( s+ T' 1 ) + e 2 ( s +T' in2 )) 



(Q - 1)(1 + e 2(l+n)(2+n 7l )) _ (Q + 1 )( e 2C(x+ 7l ) + e 2(5i+ 7l n2)) 
jNTK2/- \ _ 2(Q 2_ 1)e n(, +7l)(e 2(, +7l ^)_ 1) 

2 v. >7iJ (Q_ 1 )( 1 + e 2(i + n)(s+n 7 i))-(fi + i)(e2^+7i)+e 2 ^+7i« 2 )) ' 



This one-parametric family of solitary waves, depending on the real integration constant 71 G R, 
is the whole manifold of two-component non-topological kinks of the modeH. 

Next, Ito and Tasaki [1] addressed the issue of the kink stability of the kink solutions from a 
very interesting point of view: applying the Morse index theorem to the kink orbits they showed 
that the unique stable kink in the SI 2 < 1 range is the TK2 kink. These authors realized that 
all the other kink orbits starting from one of the vacuum points pass through one of the two foci 
(fif± = (J)i± = ±0,^2 = 0) °f the ellipse (|HJ); these points are thus respectively conjugate to one 
of the two unstable points <f>± = {4>\. = ±1, 2 = 0) because each kink orbit starting from 0_ 
(respectively 0+) passes through 4>f+ (respectively <pf~). In Reference |28j Guilarte extended this 
idea to embedding the kink stability problem in the degenerate Morse theory of the configuration 



2 In fact, there are other two parameters, 72 £ R and v G (0, 1), which are fixed by setting the kink center of mass 
system . 
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space, bearing in mind that, in this framework, the topology of an infinite dimensional manifold 
-the configuration space of the MSTB model- is determined from the critical points -the kinks- 
of a functional -the energy of the MSTB model-. A full analysis of this approach was achieved 
in Reference [29|, where the Morse's Lemma, governing the small fluctuations around the kink 
solutions, suggests how to deal with quantum corrections. This goal was fully attained from this 
point of view in Reference [30] . To finish this subsection about the history of the MSTB model, let 
us mention that in [31] the decay amplitude from the TK1 to the TK2 kink was computed in the 
quantum MSTB model by using the Witten-Smale version of Morse theory. 

Two properties of the kink variety of the MSTB model induced its early perception as singular 
and unique: a) The existence of several kinds of kink: the TK1 kink, the TK2 kink and the 
NTK2(7i) kink family, b) The relationship between its energies (the kink mass sum rule). 

1.2 The discovery of new models of the MSTB type 

The critical feature of the MSTB model allowing this perception to be overcome is the fact that its 
analogous mechanical system is an integrable system of Liouville type. In |32| Alonso, Gonzalez, 
and Guilarte proposed and studied two models of the MSTB kind with different potential energy 
densities, chosen in such a way that their analogous mechanical system would be of Type I in the 
first model and Type III Liouville integrable system in the second case. Together with a selection 
of parameters ensuring the existence of a discrete set of maxima for the mechanical potential -five 
in the first, four in the second model- this work revealed for the first time that the MSTB is by 
no means unique but, rather, a member of a broad family of models. The configuration space is 
formed by 25 (in the first case) and 16 (in the second model) topologically disconnected sectoral 
and it is expected that the kink varieties in these cases will be richer than in the MSTB system. 
This is indeed the case, and several families of stable and unstable topological kinks were discovered 
with various new kink mass sum rules arising between the different kink types. 

Model B, the second model in [32], is Liouville Type III, i.e., its analogous mechanical system 
is a two-dimensional Hamilton-Jacobi separable system in parabolic coordinates. Remarkably, 
this system arises for a particular value of the coupling constant of a model that is important in 
supersymmetric field theory. The dimensionally reduced Wess-Zumino model discussed by Bazeia 
et al. in [91 [10] is not in general of Liouville Type. Only for two special values of the coupling 
constant does the analogous mechanical system fall into this class: one of the values leads to a 
mechanical system of Type III already mentioned, and the other one gives a Liouville system of 
Type IV, Hamilton-Jacobi separable in Cartesian coordinates. The important point is that, given 
its supersymmetric origin, a solution of the Hamilton-Jacobi equation (not complete) is always 
known and trajectories associated to this solution (the superpotential) can be found. This strategy 
was used in [33J to find a whole family of topological kinks in the BNRT (acronym from |lUj . 
in analogy with the MSTB acronym from [T7] and [18]) model. The difference with model B, a 
particular case, is the HJ-separability that allows one to know a complete solution of the Hamilton- 
Jacobi equation and, therefore, knowledge of all the trajectories. The point is that for models 
with analogous mechanical system of Liouville Type, complete integrability permits an exhaustive 
knowledge of all the kink orbits, which, in turn, provides a full understanding of the solitary wave 
manifold and kink mass sum rules between them. Recently, this analysis has been performed in [33] 
for some Type II Liouville systems; the field theoretical models analyzed have analogous mechanical 
systems HJ separable in polar coordinates. 

, / N \ , 

In general, there are N + 2 ■ I j = N topologically disconnected sectors in models of this kind with TV 
zeroes of the potential energy density 
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A general pattern emerges from the study of the manifold of kink solutions in these models. 
There are two kinds of solitary waves: basic kinks with energy density localized at one point 
in R and families of composite kinks, non-linear superpositions of several basic kinks, such that 
their energy densities are localized at several points in R. The distance between the basic lumps 
of energy is determined by the parameter fixing the corresponding kink orbit. The situation is 
perfectly illustrated in Figure 1. By using a solid line in the three graphics we represent the profile 
of the three kinds of solution in the MSTB model, i.e, the TK1, the TK2 and a member of the 
NTK2 family of kinks respectively. We also depict the energy density of these solutions by using 
a dashed line. We notice that the energy density of the TK1 and TK2 kinks are localized around 
one point while the energy density for the NTK2 kink is localized around two points. The NTK2 
is composed of the other two kinks. The parameter 71 in (|lip fixes the distance between the TK1 
and TK2 kinks. In References [35} [36] a description of the evolution of this distance in time has 
been worked out for slow speeds. Although we are dealing with non-linear systems, the adiabatic 
evolution of composite kinks can be dealt with as geodesic motion in the kink parameter space with 
respect to the metric induced by the field kinetic energy density, thus providing a good description 
of kink scattering at low energies. From this point of view, kink mass sum rules become natural: 
these rules are forced by energy conservation when radiative processes are absent, i.e., the energy 
of a composite kink is the sum of the energies of the component basic kinks. 

One might wonder at this point if similar structures arise in systems with more scalar fields. In 
References |37l 139] a positive answer addressing field theoretical systems with three scalar fields was 
found. The analysis is even more difficult, because three-dimensional Hamilton-Jacobi separable 
systems not abound. Now, the composite kinks are formed by three basic kinks. Stability issues 
are also more difficult to deal with, although in [38] a thorough study of kink stability was achieved 
in the generalization of the MSTB model to three fields. The adiabatic dynamics of three slow 
moving kinks has also been recently unveiled in [40J. 

1.3 Generalized MSTB Models 

The work in the papers referred to above shows that the MSTB model is not a unique or special 
model. Several models are known to share essential characteristics with the MSTB model, in such 
a way that the kink variety of each model and its properties can be identified by a parallel analysis. 
Our goal in this paper is to study the general structure of this type of models, which we shall call 
generalized MSTB models. All these systems share the same critical properties: their analogous 
mechanical system is a Type I Liouville system having a discrete set of maxima in the mechanical 
potential energy; these MSTB-type models are distinguished by the number of zeroes -the maxima- 
of the non-positive mechanical potential energy. The kink variety can always be identified as the 
manifold of separatrix trajectories between the bounded and unbounded motion of the integrable 
mechanical system of Liouville Type I. Remarkably, the structure of the kink variety can be unveiled 
in all of these models without solving any OD equation. Moreover, the kink mass sum rules can 
also be guessed simply from the analytical form of the potential density energy. 

This paper is organized as follows: In Section §. 2 we shall define generalized MSTB models 
as field theoretical two scalar field models with analogous mechanical system of Liouville Type I. 
These two-dimensional mechanical systems are Hamilton-Jacobi separable in elliptic coordinates. 
Thus, in this Section, we shall translate the generalized MSTB models into elliptic coordinates, 
mapping the internal field space R 2 to the infinite strip E 2 = R + x I = +00) x [— O, fi], where 
0, G (0, 1) is a real number between and 1. In Section §. 3 we shall study the kink manifold 
that arises in this kind of model. In order to do this, the distribution of the zeroes of the potential 
energy density is crucial. In generalized MSTB models these zeroes are located at the nodes of a 
grid or reticulum in the "elliptic strip" E 2 . The distribution of zeroes of the potential determines 
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the classification of generalized MSTB models in four types. Generically, the reticulum is formed 
by gluing several cells where the kinks are confined. The zeroes of the potential energy density 
live at the vertices of the cells and/or at the foci of the elliptic and hyperbolic curves that set the 
curvilinear coordinate system. There are three types of cells distinguished by the nature of the 
vertices. The type of cell determines the behavior of the kink family confined in the cell. In general, 
there are basic and isolated kinks living in the bars of the grid surrounding the cell and families 
of kinks living inside. Each member of one of these families is a non-linear superposition of basic 
kinks: these are composite kinks and have their energy density localized around several points, 
each of these component lumps having the same energy density as a basic kink. In this Section, we 
shall also explore new and generalized kink mass sum rules and we shall investigate the relationship 
between the integrability of the analogous mechanical systems and the existence of superpotentials. 
In Section §. 4 we shall describe the kink variety of a specific model. The MSTB model [T7] involves 
a quartic algebraic potential energy density, whereas the potential energy density of model A in 
is a polynomial in the fields of sixth degree. In order to deal with completely new generalized MSTB 
models, we shall analyze a model characterized by a eighth degree polynomial in the two fields as 
the potential energy density. The polynomial depends on two coupling constants in such a way 
that a process of bifurcation arises: the number of elements in Ai, static homogeneous solutions of 
the field theoretical model or the fixed points of the Newton equations of the analogous mechanical 
system, depends on the values of the coupling constants. There are three regimes where we shall 
study the kinks and also the kink mass sum rules. Finally, the summary of the main results and 
the outlook on future research in similar models are offered in Section S. 5. 



2 The structure of generalized MSTB models 
2.1 The configuration space 

We shall study field theoretical models of two real scalar fields defined in a (1+1) Minkowskian 
space-time, whose dynamics is governed by the action functional: 



sy>] 



dxdt 



1 2 



0=1 



(12) 



We use Einstein sum convention for Greek scripts and choose the metric tensor to be: goo = —gu 
1) 5oi = 5io = 0- By eft = X^a=i $a£a-, a = 1, 2, we denote the scalar fields: 



<f>{x Q ,xi) = <f>i(x ,xi)ei + fa(x , xi)e 2 : 



jl.i 



The vectors e form an orthonormal basis in R 2 and U(fa, fa) is the (non-negative) potential energy 
density, which vanishes in a discrete set of points, M.. For the sake of simplicity we assume that 
the magnitudes involved in (|12|1 are dimensionless. 

A point in the configuration space of the system is a configuration of the field of finite energy; 
i.e., a picture of the field at a fixed time such that the energy E, the integral over the line of the 
energy density, 



E[4>] 



dx£[4>(x)] 



1 



dx 



dx 



(13) 



is finite. Therefore, the configuration space is the set of continuous maps from R to R 2 of finite 
energy: 

C = {<f(x) E Maps(R,R 2 ) / E[<f>] < 00} 
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In order to belong to C, each configuration must comply with the asymptotic conditions 



lim 4>(x) G A4 , lim -f- = (14) 

x— >±oo x— >±oo ax 

where .M is the set of zeroes of the potential energy density U((j>), 

M = G M 2 C Maps(R,M 2 )/[/[0 (/) ] = , I = 1, 2, ■ ■ • ,N} 

We shall also choose U(4>) non- negative and such that Ai is also the set of absolute minima of 
U, i.e., the set of the static homogeneous solutions of the model because the Euler-Lagrange field 
equations of (fT2"j) 

d 2 6 a d 2 6 a dU , . 

are satisfied by the elements of Ai. Moreover, bearing in mind that time evolution is continuous -a 
homotopy transformation-, if Ai is a discrete set the configuration space is the union of disconnected 
topologically sectors C IJ : 

C= |J C IJ . 

I,J£M 

The configurations in C IJ are characterized as those which respectively tend to the elements 
and (f>( J ^ of Ai at minus and plus spatial infinity. 

According to Rajaraman [6]: A kink is a non-singular solution of non-linear field equations M5\) 
whose energy density, as well as being localized, has space-time dependence of the form: e(t, x) = 
e(x — vt) ; where v is some velocity vector. Because we are dealing with a Lorentz invariant system, 
solutions with a temporal dependence like this are obtained from static solutions by means of a 
Lorentz velocity transformation: L(v)(f)(x) = <p L (t,x) = (j) ^-^=4^. Bearing this in mind the PDE 
(fl~5|) can be restricted to the ODE 

d^a = dU_ a= _ fl = y U (16) 

dx 2 d(j) a ' ' dx 2 

in such a way that the search for some kinks or solitary waves is tantamount to studying an 
analogous mechanical system, where we interpret the static scalar field 4>(x) = ^i(a;)ei +(j)2(x)e2 as 
the particle position vector; x as the particle time, and U = —V as the particle potential. Moreover, 
the action and the Lagrangian of the mechanical system respectively play the roles of the energy 
functional and the energy density in the field theory, see (|13p . In sum, solitary waves or kinks 
are derived from stationary points of the functional energy (|13|) that comply with the asymptotic 
conditions (|14p . Therefore, kinks asymptotically connect two elements of Ai or, homogeneous 
solutions of (fT6l). 



2.2 Algebraic Type I Liouville potentials 

The potential energy density of the MSTB model is a positive quartic algebraic expression in the 
fields, see (0). In this sub-Section we shall characterize all the algebraic (polynomial) potential 
energy densities that have Type I Liouville systems as the analogous mechanical model. Because 
the polynomials will in general be of higher order than four, the cardinal of Ai may be greater than 
in the MSTB model, promising richer kink varieties. 
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Proposition 1: All the members in the (q + 2) ^-parametric family of models with algebraic 
energy density of the form: 



n=0 j=l 



2Vt(p\ cos 



J7T 

n+ 1 



(17) 



a polynomial in the (pi and (p 2 fields of degree 2q, have analogous mechanical systems that are of 
Liouville Type I, RTF, and, therefore, these mechanical systems are Hamilton- J acobi separable in a 
system of elliptic coordinates. 

The proposition will be proved later in this sub-Section. In formula (|17p b n , n = 1, . . . ,q and 
SI are the parameters or coupling constants distinguishing between the members of the family. 
In order to support kinks /solitary waves, or non-trivial separatrix trajectories in the analogous 
mechanical system, the choice of b n parameters must be restricted in such a way that U((pi,(p2) 
must be bounded below by a discrete degenerate set of constant field configurations that are the 
absolute minima of U ((pi, 02 )• Also, some transformations, e.g. scaling of the space-time = Xx^, 
scaling and rotations in internal space, (p a = J(j> a , (pi = cosa0i — sina02, 02 = sina0i + cosa02 
lead to equivalent models, allowing us to restrict the values of three parameters. 

Arguing as Magyari and Thomas for the MSTB model in [25], we notice that the analogous 
mechanical system for the potential energy density (|17p is completely integrable: the static field 
equations (|16p admit two independent invariants in involution: 



1 2 { dx J 2 \ dx , 

71=0 J = l 



dx 



dx 



dx 



-2fi<£i cos — + n 2 J + n 2 Yl ( <P\ + 



'1 \ 2 o2 ( d< p2^ 1 



n=0 



2S701 COS + O' 

n + 1 



71-1 

*1 11(^1 + ^ 

i=i 



i>2 — 2001 cos 



J 71 " 
n + 1 



+ S7< 



(18) 



(19) 



The first invariant is the energy of the mechanical system and the second one is also quadratic 
in the velocities, a certain deformation of a first-integral of Runge-Lenz type. Therefore, the field 
theoretical model is associated to an analogous mechanical system of Stackel type. By choosing 
an appropriate system of coordinates dynamical systems of this kind are always Hamilton-Jacobi 
separable. The form of the second invariant dictates that elliptic coordinates will produce a sep- 
arable Hamilton-Jacobi equation (Type I Liouville systems) and all the solutions can be found by 
quadratures. In the MSTB model, this property was used by Ito [27J) to find all the solitary wave 
solutions. 

To study algebraic Type I Liouville systems we define Euler elliptic coordinates by means of 
the map: 

^ : E 2 = x [SI, 00) — > M. 2 /Z 2 

(v,u) — > (01, 2 ) 



&(«) 



yj {(p\ + SI) 2 + <p\ — y/ {<pi — SI) 2 + (j>; 



2 ' >i2V ' 2 

which is a diffeomorphism between the elliptic open strip (— S1,S7) x (S7,+oo) and the Cartesian 
upper half-plane 02 > 0. Explicitly, the inverse map from M 2 /Z 2 to E 2 , is: 



h uv 



,) = ±lvV-S7 2 )(S7 2 -*; 2 ) 



(20) 
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such that the coordinating curves are the hyperbolae and ellipses 



r , ~r 



V 2 O 2 — V 2 U 2 U 2 — O 2 

whose foci are the points: F± = (±0, 0). It is necessary to use two copies of the elliptic strip E 2 
to cover the whole Cartesian internal plane: £^ maps E 2 into the fa > half-plane whereas £^ 
maps the other copy of E 2 into the (j)2 < half-plane. Translation of information from E 2 to the 
Cartesian R 2 internal plane must take this fact into account by smoothly gluing data across the 
border between the two copies. In particular, the straight lines u = O and v = ±0 in E 2 determine 
the axis (j>\ in R 2 ; u = O is the (j>\ £ [—0,0], (f>2 = interval and v = ±0 are respectively the 
<f>\ £ [— oo, —0] (- sign) and [0,oo],(/>2 = (+ sign) half-lines, all together forming the abscissa axis 
in M?. Note that the parameter setting the foci of the coordinating curves is chosen to be the 
coupling constant O in (fT7|) . 

Applying the map (f20|) to the potential energy density (fT7|) we obtain: 



ftl7(&,fc) = -^-^j^b n {u 2n+2 -v 2n+2 ) = ^-^(f{u)+g(v)) (21) 

n=0 

/:[fi,oo)->R ; u m- /(u) = £ & n u 2n+2 (22) 

?1=0 

5 : [-0, 0] -> R ; v ^ = - M 2n+2 . (23) 

n=0 

In elliptic coordinates the potential energy density (|2ip is the product of the u ^_ v i factor times the 
sum of two polynomials, which depend only on one of the elliptic variables, -gj^-g is in all places 
positive in the strip E 2 except at the foci F± = (u,v) = (±0,0) (the points (<j>\,(j>2) = (±0,0) 
in R 2 ), where the denominator vanishes. Despite this, the potential energy density is not singular 
at the foci F because the tendency to zero is at least as strong in the numerator of (|21[) as in the 
denominator. This can be seen in the factorization: 

u 2n+2_ y 2n+2 _ ( U 2 - V 2 )( U 2 -2 U V COS — \~V 2 ) (ii 2 — 1uV COS — — — hf 2 ) •. . . ■ (it 2 — IllV COS — — — \-V 2 ) . 

y n n + 1 n + 1 ' v n + 1 ' 

(24) 

which also explains the symmetric form of the functions / and g, i.e., g(z) = —f{z). The inverse 
image in Cartesian coordinates of (|2ip provides the strange form of the potential energy density 

m u q {4> x ,(j>2). 

The main point is that the static energy functional (|13p in elliptic variables reads: 



fi*l?U ^ 1 I A l u2 ~v 2 f du Y , l u 2 -v 2 (dv\ 2 f(u)+g(v) 

The functional £^.E7[<^i, <fo] : Maps(R,E 2 ) — > R can be interpreted as the functional action for a 
particle moving in E 2 under a mechanical potential energy V[u(x), v(x)] = — £qU[4>i(x), 4>2(x)] if 
the particle time is x. In this conceptual framework, the change to elliptic coordinates induces a 
flat, but non-Euclidean, metric: 

ii 2 — ii 2 ii 2 — ?; 2 

9 UU M = ^—- , g m (u,v) = - 2 j , g^( u ,v) = = g" u (u,v) . 

u z — \ l z \l z — v z 
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In any case, it is well known that mechanical systems with dynamics ruled by an action such as 
(j'25|) are Type I Liouville systems, see [41j . Therefore, there are two first integrals: 



t* t _ 1 u 
- ^ 



2u 2 -n 2 \dx J 



du\ 



1 u 2 



u 2 ) 



n 2 -v 2 fdu\ 



— 



u 



20 2 
o 2 



*A _ f(u) +g(v) 
dx J u 2 — v 2 



(26) 



n 2 \dxj o 2 



dx J 



(u 2 -n 2 )g(v) | (O 2 - t; 2 )/ (^ 



U — ir 



which allow us to solve the motion equations. Back in Cartesian coordinates, (|26p and (|27p become 
respectively (fl"8|) and (fT9|) . fully proving the Proposition. 



2.3 Generalized MSTB models 

Generalized MSTB models are a subset of those models having algebraic Type I Liouville systems 
as analogous mechanical system. In order to characterize with precision this kind of model we 
recall how the MSTB model arises between the q = 2 algebraic Type I Liouville systems: 

U q=2 (<l>i,<h) = b + h (V? + 0|-2O0i cos | + 2 

(it \ ( 2tt 

0? + 4>l - 2O0i cos - + O 2 j ( <j){ + <j) 2 2 - 2O0i cos — + O 2 

= ai(0i + 02 + a 2 ) 2 + a30i + «40| 

The criterion is that, up to equivalences, U q= 2 must be non-negative, having a discrete set of 
zeroes with more than one element; these properties give a hint of the existence of solitary waves. 
Demanding U q= 2 > with a discrete set of minima requires a\ > 0, > and 04 > 0. Also a\ can 
be fixed by re-scaling the space-time, a>2 by re-scalings of the fields and 03 by rotations of the fields. 
In the MSTB model, we choose a% = 5, 02 = — 1, a 3 = such that only a positive parameter, 04, 
is left and there are two zeroes of U q= 2- 

Another model belonging to this kinship has been dealt with in the Literature. The q = 3 
family of potential energy densities is: 

U q =3 (01, 02) = b + 61 (0? + 4> 2 2 - 2O0i cos I + O s 

+ ^2 (0? + 02 - 2001 COS I + O 2 ) ^0 2 + 02 - 2001 COS ^ + ^ 

+ 63 (0 2 + 02 - 2001 COS j + O 2 ) (^t)\ + 02 - 2001 cos ^ + 2 ^J ^0? + 02 - 2001 cos ^ + O 2 

= ai(0 2 + 2 ) 2 (0? + 2 + a 2 ) + (a 3 2 + a 4 )(0 2 + 0l) + ^0102 + 03)^ + 

4ai 

?7 9= 3 is a five-parametric algebraic polynomial in 0i and 02 of sixth order that can be understood 
as a deformation of the Chern-Simons-Higgs potential energy density: 

C/CSH(01,0 2 ) = (0? + 0|)(0? + 0! " I)' 

The A model investigated in Reference is a special member of this family of q = 3 with five 
zeroes of U q= ^. Therefore, Ai has five elements in the A model and the configuration space is the 
union of 25 disconnected topological sectors. 
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The potential energy densities in the MSTB Q and A model [36] in elliptic coordinates take 
the form of formula (|2ip by choosing f(u) and g(y) respectively to be: 



/(«) = \{u 2 ~ 1) V - O 2 ) , g(v) = \{v 2 - 1) 2 (0 2 - v 2 ) 

f{u) = \u 2 {u 2 -l) 2 {u 2 -tf) , g(v) = 1 -v 2 (v 2 -l) 2 (n 2 -v 2 ) . 

Definition 1: Generalized MSTB models are (1+1) -dimensional field theoretical models of two 
scalar fields such that their potential energy densities have the following form in elliptic coordinates: 

n+m n+m 

f(u) = (u 2 -n 2 )u 2a ° H (u 2 -a 2 ) 2a > , g(v) = (n 2 -v 2 )v 2a ° J] {v 2 -a 2 ) 2a > . (28) 
j=i 3=1 

Here ao,ctj,n,m 6 N are natural numbers and there are n + m + 1 coupling constants: o~j,j = 
1, 2, ■ • • , n + m, and 0. With no loss of generality, we order the coupling constants in such a way 
that: 

< o\ < a2 < ■ ■ ■ < cr n -i < a n = 0, < a n+ \ < ■ ■ ■ < cr n+m 
This choice is guided by the following criteria: 

• The potential energy density is positive for all the values of the coupling constants because 
all the terms involve even-power factors. Notice that the factors {u 2 — Q 2 ) and (J7 2 — v 2 ) 
in the functions f{u) and g(v) are semi-definite positive because of the range of the elliptic 
coordinates u G [O, oo] and v G [—0,0]. The presence of these factors allows us to comply 
with the condition f{z) = —g(z), preserving the positiveness of the potential energy density. 

• The functions f(u) and g(v) have been also chosen in order to maximize the cardinal of the 
set of minima or zeroes Ai. For this reason these functions involve only real roots. In this 
way the models will have a richer kink variety . 

Because solitary waves come from separatrix trajectories complying with the asymptotic conditions 
(|14p . kinks will arise when = = ^^2- Indeed the condition = (|26p is tantamount to 
the virial theorem and the kink energy density reads: 

2 r „. . . u 2 — v 2 ( du\ 2 v? — v 2 f dv\ 2 , . 

«fc "1 = jfz? [/<») + »M1 = ^ttp (s j + W^P U ) (29) 



2.4 The structure of the set M. of zeroes of [/(</>i,0 2 ) 

Bearing in mind that solitary wave solutions asymptotically connect two homogeneous minima of 
U (4>\ ,<p2)i m this sub-Section we shall study the distribution of elements of Ai in the internal plane 
M 2 . We shall show that the minima are located at the nodes of a rectangular mesh or reticulum 
both on the elliptic strip and on the Cartesian plane. This grid comprises a number of cells that 
depends on each specific model. Kinks are confined inside the cells and their qualitative behavior 
is determined by knowing only the type of cells where they are constrained. Moreover, on the 
edges of the cells single kinks emerge, while inside the cells families of composite kinks (non-linear 
combinations of single kinks) arise. 

In order to systematize the analysis of generalized MSTB models we shall start by describing 
the set M. of zeroes of £^£/(0i, fa) on the elliptic strip E 2 . From (|2"Tj) . one notices that the zeroes of 
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the potential energy density are roots of the functions f(u) and g(v). We shall denote respectively 
by 

M g = {v G [~n, Q] I g(y) = 0} and Mj = {u£ [fi, oo) / /(«) = 0} 

the sets of real roots of ^(u) and /(«)• We stress the following points about the elements of M. g 
and M /: 

• If dj is a root of 5(f), — <Jj is another root, see formula ([28]) , We shall simply denote: 
a -i = -°y 

• Although f(u) and g(u) are almost identical functions their domains are different: the domain 
of g is a finite interval, Doing = [— fi,fi], whereas the domain of / is the infinite half-line, 
Dom/ = [0,oo). For this reason, in general the numbers of roots of these functions are also 
different. We shall denote by n and m + 1, respectively the number of positive roots of g(v) 
and f(u). Therefore, card.M/ = m+1 and, either card7W 9 = 2n + l if v = is a root of g(v), 
or, cax&Mg = 2n if is not a root. We shall show later that the pair (n,m) characterizes 
the qualitative behavior of the kink variety in the generalized MSTB models; models with 
the same pair (n, m) have kink varieties of the same type. 

• It is convenient to order together the roots crj of the two functions g(v) and f(u) from the 
lowest to the highest. Let J be the set of labels: J = {— n, — (n — l),...,n — l,n}. The 
set of roots of g(v) is: M. g = {aj G [— j G J}. If g(v) has no zero root j = ^ J. 
Analogously, let if be another set of labels: K = {n, n + 1, ...,n + m}. The set of roots of 
f(u) is: Mf = {&k G [il,oo) / A; G K}. It is convenient to use the global set of labels J U K 
because j G JUK = {— n, — (n — 1), ...,n,n + 1, + m} jointly label the roots of the two 
functions g(v) and f(u) (see Figure 2). 

• It is clear that v = ±f2 are the lowest and highest roots of g{v) whereas u = VL is the lowest 
root of f(u). Therefore, the (increasing) ordering of the positive or zero roots of M g and A4f 
is as follows: 

< o\ < cr 2 < ... < cr n _i < cr n = Q < a n+1 < ... < cr n+m 



-4 Dom g( v) >■ *Dom f(u) + 

0-„=-Q. o_i o = 0Oi a„=Q. a n+ i a n+ ,„ 

I 1 1 1 1 1 1 

i =—n -1 0? 1 n n+1 n+m 

Figure 2: Ordering of the roots of the functions g(v) and f{u). 



The set Ai of zeroes of U(cj>i, <f>2) has a straightforward connection with the zeroes of the sum of 
f(u) and g(v), (the numerator in (f2"Tj) ). in E 2 : 

M f+g = M g x M f = (J (v = a j ,u = a k ) 

or, back in M?: 

M f+ g = &M f+g = |J U =%p;<h = ±y(o* - - cr|) 
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The elements of M. are obviously zeroes of /(u) + g(v) but the converse statement is false: M. C 
Aif +g . The reason for this is that the foci F± of the coordinating curves belong to -Mf+ g but the 
potential energy density (f2Tj) might not vanish at these points because of the 2 factor in ([2~T]) . 

2.5 Types of generalized MSTB models 

There are two pairs of excluding possibilities, (1) or (2), (a) or (b), prompting a classification of 
the generalized MSTB models: 

- (1) v = ±0 are simple roots of the function g(v), or, (2) v = ±0 are multiple roots of the 
function g(v). 

- (a) The function g(v) has no zero roots, or, (b) the function g{v) has zero roots. 
Therefore, there are four kinds of generalized MSTB models: 

• TYPE I-la: All the roots of g(v) are non zero and v = a± n = dbfi are simple roots. In this 
case the foci F± do not belong to the set £q*M = M of minima of the potential energy 
density in E 2 and: 

M = M f+g -{(±n,n)}= |J {(a j ,a k )}-{(±n,n)}= [j v j}k . 

(j,fc)eJxK (j,fc)eJxK 

Note that J = {— n, —1, 1, ... , n}, i.e., the element j = is not accounted for in the set of 
labels J because there is no zero root of g(v). The number of elements of the minima of the 
potential energy density is straightforward in E 2 : card(A4) = 2n(m + 1) — 2. To give the 
corresponding number in M 2 is a tricky business. In general, to each minima in E 2 correspond 
two minima in R 2 : 

£(<r±j,vk) = (^p,±^(a 2 - tt 2 )(^ 2 - a 2 ) ) = v% , k = n + 1, n + 2, • • • , n + m 

This is only true in the interior of E 2 . Each minima living at the boundary maps only in 
one minima in 1R 2 sitting on the abscissa axis: ^Q(a±j,Q) = (a±j,0), j = 1,2, ••• ,n — 1, 
£q (±f2, Ufc) = (±<Tfc, 0), k = n + 1, n + 2, • • • , n + m. Therefore, the number of minima of the 
potential energy density in R 2 is: card(A^) = 4nm + 2(n — m) — 2. In Figure 3 the distribution 
of all these points is shown, both in E 2 and M 2 . The straight lines joining the minima in the 
elliptic strip (hyperbolae and ellipses in the Cartesian plane) form the reticulum associated 
with this type of model. 

u 

• • • • • • # • a„ + , 



1 p — 1 1 — 1 1 1 1 1 1 1 1 — 1 1 — 1 1 c n+2 

II (i 1 1 < 1 (i i 1 < t a„ + 1 

F (p n I n (| M fI V 

Figure 3: M (points) and the reticulumlZ (solid lines) in E 2 (left), M. (points) and the reticulumlZ (solid 
lines) in M 2 (right) for the Type I-la models. 
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• TYPE I-lb: In this case v = o~o = is a root and a± n = ±J1 are simple roots of g(v). The 
foci F± do not belong to M. in this type of model either and: 

M = M f+g -{(±n,n)}= |J {(^,<T fc )}-{(±^)}= |J « j)fc , 

(j,fc)eJxK (j,fc)eJxK 

where now J = {— n, —1, 0, 1, n} because v = is a root of g(v). The number of elements 
of M. is: card(TW) = 2n(n + 1) + m — 1. In the count of minima in M 2 we must add to the 

minima of Type I-la models the points: £q (0, at) = (0, iy/fr 2 — il 2 ), k = n, n+1, ■ ■ ■ ,n+m, 
all of them living on the ordinate <f>i = axis. Therefore, card(.M) = Anm + 2n — 1 because 
all the minima in the v = line of E 2 go to two minima in M. 2 , except the (k = n,j = 0) 
minimum which is mapped at the origin ofM 2 . InFi gure 4 the minima of the potential energy 
density in this type of model is depicted both in the elliptic strip and in the Cartesian plane. 
The reticulum of straight lines joining these points is also plotted. Different cells with respect 
to the Type I-la cells arise close to the ordinate axis due to the new Type I-lb minima. . 




Figure 4: A4 (points) and the reticulum 1Z (solid lines) (left) , M. (points) and the reticulum 1Z (solid lines) 
(right) for the Type I-lb models. 

• TYPE I-2a: In this type of model v = a± n = ±Q are multiple roots of the function g(v), 
such that the foci F± belong to the set M whereas v = is not a root. Therefore, 

M = M f+g = \J {(a j ,a k )}= \J v jtk 
(j,fc)eJxK (j,fe)GJxK 

where J = {— n, —1, 1, n}. The number of elements in M is: card(TW) = 2n(m + 1). In 
M. 2 , however, this number is: card(A^) = Anm + 2(n — m), because the two new minima at 
the foci in E 2 are one-to-one mapped to the foci in M 2 , see Figure 5. 




Figure 5: M (points) and the reticulum 1Z (solid lines) (left) A4 (points) and the reticulum 1Z (solid 
lines) (right) for the Type I-2a models. 

• TYPE I-2b: In the last type, v = o"o = is a root and v = a± n = ±$7 are multiple roots of 
the function g(v). There are minima of the potential energy density in E 2 at the foci F± and 
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on the v = line. Therefore, 



M = M f+g 



U {(Oj,<J k )} 

(i,fc)eJxK 



U Vj, k 

(i,fc)eJxK 



where J = {— n, — 1, 0, 1, n} and: card(TW) = (2n + l)(m + 1). Passing to Cartesian 
coordinates, each minimum at the two foci and the (v = 0, u = CI) point go to one minimum 
in the <j>2 = abscissa axis, but the m remaining new minima (v = 0,u = a/~), k = n + 
l,n + 2,---,n + m, have two image minimum in M 2 . Accordingly, the count of minima in the 
Cartesian plane is: card(.M) = 4mm + 2n + 1, see Figure 6. 



#—4 



<J„ + i 



CT-l 



CJl 



<J„-i C5„ 




Figure 6: M. (points) and the reticulum 1Z (solid lines) (left) , M (points) and the reticulum 1Z (solid 
lines) (right) for the Type I- 2b models. 



In all four cases, the network in E 2 is formed by the horizontal segments 

= {{v, u) e E/crj < v < <Tj+l,U = CTfc} 
joining the points (<7j, a^), and (cTj+i,^) and the vertical segments 

r j'j +1 = n ) G E l v = a ji a k < u < o- fc+ i} 

which join the points (aj,ak) and (<jj, crfe + i). Recall that j G J, A; G Kandj + 1 means the following 
element in the set J. The reticula 1Z in R 2 become the sets of orthogonal hyperbolae and ellipses: 




i G JUK 



2.6 Types of cells 

The reticulum 7?. on the elliptic plane is formed by rectangular cells, see Figures 3-6. We shall 
denote by C^'j^ the cell whose vertices are the elements of -Mf+ g ((Tj,o~k), (o~j,ak+i), (aj+i,o~k) 
and (cTj + i, (Tk+i)- Therefore, the edges of this cell are the segments r^ fc +1 , fjf +1 , ^j+i^+i and 
r j~j+i +1 - We distinguish the following types of cells: 

• Cells of Type 1: The four vertices of the cell are elements of A4, see Figure 7a. This type of 
cell arises in all the kinds of generalized MSTB models and it is the usual pattern for sufficiently 
high values of m and n, see Figures 3, 4, 5 and 6. 

• Cells of Type 2: Three of the vertices of the cell are elements of M whereas the fourth vertex 
is one of the foci of the model, see Figure 7b. The cells C^'^^ ±n that arise in the type I-la and 
I-lb models are of this type. In this case the foci F± are not zeroes of the potential energy density. 
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• Cells of Type 3: Only two of the vertices of the cell belong to M whereas the other two are 

1 2 

the foci of the coordinate system, see Figure 7c. This type of Cj 1 1 cell arises in the special case 
(n = l,m) of Type I-la models. Note that (n = l,m = 1) is precisely the MSTB model. 
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Figure 7: Types of cells: (a) Type 1 cell (the four vertices are zeroes of^U), (b) Type 2 cell (three of the 
four vertices are zeroes of £ n U ) and (c) Type 3 cell (two of the four vertices are zeroes of S,qU). The edges 
of the cells have also been plotted as well as the singular Kinks living on them. 



In sum, Type I-2a and I-2b models exhibit only cells of Type 1, Type I-la and I-lb models have 
cells of both Type 1 and Type 2. Cells of Type 3 arise only in (n = 1, m) Type I-la models. 



3 The structure of the kink varieties 



This Section is devoted to analyzing the general structure of kink varieties in the generalized MSTB 
models. In order to alleviate the notation we shall denote the elements of M in a generic cell as 
A = B = fj+i,fe+i, C = Vj^ and D = vj + i t k, see Figure 7. The notation for solitary waves 

will be KJj where the superscripts designate the elements of M (or M) that are asymptotically 
connected by the kink solutions, and N stands for the number of energy lumps supported by 
the traveling wave. Bearing in mind the invariance of generalized MSTB models under spatial 
reflections, x — > — x, one may identify Kjtf and Kj!^'" as a single type of solution. Also, we shall use 
the term single kinks to refer to those solutions whose energy density is localized around a single 
point (N = 1) whereas we use composite kinks to refer to those solutions whose energy density can 
be localized at several points (J\f > 1). 

There are two main categories or kinds of solitary wave solutions: (1) Singular (isolated) kinks, 
which are the single or basic kinks of generalized MSTB models. (2) One-parametric families of 
composite kinks, which are non-linear combinations of the single kinks. Composite kinks are formed 
by several lumps whose separation is determined by the value of the parameter of the family. 



3.1 Singular kinks 

Proposition 2. There exist kinks whose orbits are located on the border of the cells Cj'^i, 
j 6 J, k G K ; i.e., on the segments or r k ^ +1 of the elliptic strip E 2 (pieces of hyperbolae 

and ellipses in the M 2 Cartesian plane). For reasons to become clear later we shall refer to these 
solutions as singular kinks. 

Proof: The critical value = = Chh, where the two invariants cease to be independent, is 
reached by all the trajectories or orbits satisfying the system of first-order ODE's: 



du_ V2(u2 - Q2 )/(m) dv _ ^/ 2 (W-v*)g(v 

dx u 2 — v 2 dx v? — v 2 
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see ([26]) and ([57]) . We try the orbits u = ak, k £ K (the segments or t> = Oj, j G J (the 



segments fjj k+1 ) into the (|30|) equations. Alternatively, one of the two equations is satisfied if either 
*i = or £ = 0, and the other may be immediately integrated by means of a quadrature. In sum, 
there are 

— Kink orbits living on pieces of hyperbolae r j'j +1 '- 



v(x) = cjj , u(x) = h 1 (x) , h(u) = I — 3 ' (31) 



III 



2 



(Jj)du 



Kink orbits living on pieces of ellipses r^'*, , : 



v(x)=h~ 1 (x) , u(x) = (Tk , h(v) = / — k 



(32) 



and the Proposition is proved. 

The type of singular kink, however, depends on the type of cell where singular kink orbits arise. 

• In type 1 cells there are four kinds of kinks. The four vertices A, B, C and D in these cells 
are elements of M.: A = Vj t k+i, B = fj+i^+i, C = Vj^ and D = see Figure 7a. Kink 

orbits connect two adjacent corners through the edges of the cell. Thus, singular kink orbits 
are the segments (u = &k+i) and (u = a^) respectively joining A to B and C 

to D. In the M. 2 Cartesian plane these kink orbits become elliptic arcs. 

There are also singular kink orbits connecting A to C and B to D along the segments ^j C+1 
(v = Uj) and r^ 1 + j L +1 (v = (arcs of hyperbolae in the M 2 ). We shall denote these 

four kinds of kinks as K^ B ■, K^ D , K± and Kf D . Their associated profiles are respectively 
obtained from the quadratures (f3Tj) and (f32|) . 

The main feature of this type of kink is that they are single kinks. To show this, note that: 

— The kink solutions are monotonic (increasing or decreasing) functions u(x) and v{x) of 
x. In the first-order equations (|3(J|) one immediately realizes this point because 7^ 
and 7^ for every point in E 2 except at the zeroes of f(u) and g(v), which are reached 
only asymptotically (at x = ±oo). 



The kink energy density (|29p becomes a polynomial in a single elliptic variable for this 
type of kink orbits. For instance, for v = Oj orbits the energy density reads: 

£[u,<Tj] = 2u 2a «(u 2 - a 2 ) 2 " 1 • . . . • (u 2 - a 2 _ x ) 2a ^{u 2 - a 2 ) 2a >-\u 2 - a 2 +1 ) 2a ^ 

• (u 2 - al^-^u 2 - n 2 )(u 2 - a 2 n+1 ) 2a ^ • . . . • (n 2 - a 2 n+ J 2a "+™ . 

Because this polynomial p(u) = £[u,cTj] is a semi-definite positive function in the u £ 
Wk^k+i] range that vanishes only at u = and u = (Tk+i, Rolle's theorem ensures 
that p(u) reaches a maximum. To show that this maximum is the only critical point in 
this interval requires a little more work. In fact, it is better to consider the polynomial 
p(u) as a function of z = u 2 : 

p(z) = n£+ m (z - Zi ) n > , z = , Zl = a 2 , i = l,2,---,n + m 
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The multiplicities rtj of the roots Z{ ar^f|: n{ = 2a>i except no = ao, rij = 2a j — 1, n n = 1. 
Let us consider the log p(z) function and its derivative: 

n+m /, \ n+m 

rj(z) = logp(z) = ^ n i l °g( z ~ Zi) , v'i z ) = —FT = ^2 

i=i p[z) i=o z Zi 

Excluding the points zi, z%, ■ ■ ■ , z n +m of the half-line [1, +oo) the zeroes of rf and p' are 
the same. We consider a variation of z asymptotically connecting a 2 with cr| +1 . Because 
rf' {z) = — X^i=o™ (z-zrf * s negative, r]'(z) is a decreasing function of z. This, together 
with the limiting values 

lim rj'(z) = +oo , lim rf{z) = — oo 

means that rf intersects the abscissa axis at only one point and, consequently, p(z) has 
only one critical point. Therefore, the kink energy density is localized around one point 
and this single kink is a basic lump or particle. 

In Type 2 cells, C^J^}^ _j_ n three of the vertices are elements of A4 but the fourth vertex is 
one of the foci F±; A = i>±( n _i) jn 4-i, B = v± ri . n+ \, C = w±( n -i),n an d D = F±. Kink orbits 
K± B and K^ c respectively joining A to B and A to C through the segments u = <r n +i and 
v = a n -i, see Figure 7b, are identical to those in Type 1 cells, providing similar basic kinks. 
Additionally, there is a new kind of kink orbit connecting B and C along two edges of the cell: 
the orbit starts from C, follows a segment where u = O, crosses one of the foci F±, runs along 
a v = ±Q segment, and, eventually, arrives at B. The path is the union of the trajectories 
r 



n— l.ri 



U rn,n or _ ra U r™'™_ n , an interval on the <p\ axis in the M? Cartesian plane, 

such that the solution in M? is continuous and differentiable. 

Because the orbits travel two edges of the cell, one guesses that this new type of kink will 
be combinations of two basic kinks. Therefore, we shall term them FC . Let x c be the 

point on the line where the kink profile crosses the focus. Recall that in this cell the range of 
elliptic variables is u G [f2, cr n +i] and v = [a n -\,£H\. Thus, the energy density reads: 



£[K^ FC ](u,v) 



n+m 



\\ {v{xf - a]) 2 *** if x<x c 



j=l,jytn 
n+m 



] ] (u(x) 2 - a]) 2 "* if x>x c 



Taking into account that the polynomial p(z) = n™j 1 ™_^ n (z 2 — cr 2 ) 2ai has roots at the values 
z = <J n -i and z = cr n +\, we expect a relative maximum in the interval [oW- 1> °n+l]> which is 
achieved either in the interval [cTn-i, fi] or in [f2, <T n +i]. An identical argument to that applied 
to the other type of kink shows that the maximum is the only critical point of the energy 
density. Thus, this kind of solution exhibits one lump, which may seem contradictory to the 
above guess about the composite nature of these kinks. We shall confirm later that these 
solutions are formed by two overlapping basic lumps by observing the kink mass sum rules 
in the following Section. 

In Type 3 cells of (n = 1, m) Type I-la generalized MSTB models only the vertices A = u_i 2 
and B = v\2 are elements of A4, whereas the other two vertices are the foci F± of the model. 



4 Strictly, zq — is not a root because z £ [1, oo), but this does not affect our argument. 
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There are two kinds of kinks: the kinks akin to the kinks in Type 1 cells running on 

the AB edge, the segment r^j, and the j^bf+M f + F-\j F_A j^^g The orbit of this new 
kind of kink is formed by the other three edges of the cell, see Figure 7c. The kink orbit 
starts at the point A, following the straight line v = — fi. After crossing the focus F- the 
orbit runs through the u = segment and reaches the focus F + . Finally, the kink orbit 
arrives at the point B, traveling the segment v = O. In sum, this three-stage kink orbit is: 
r 2 '| _jU ri'| 1 U rj' 2 . In the Cartesian IR 2 plane, this trajectory is the [A, B] interval of the <j)\ 
axis. Denoting by x± the points in the spatial line where the kink orbit in the elliptic strip 
E 2 passes through the foci F± we write the energy density of kinks of this type in the form: 

if x < x_ 
if x_ < x < x + 
if x + < x 

One can argue that this energy density reaches a maximum value at only one point inside the 
orbit whereas it tends to zero at both endpoints along the lines shown for the other singular 
kink orbits. In this bizarre case the kink is not a composite of another two kinks because the 
end points are shared with the only other kink in the cell. Thus, despite being built from 
three stage orbits this kink is truly a basic kink. 

In this Section we have shown that there are kink orbits living on the border of each cell C^'^j 1 . 
It remains to decide whether or not there exist solitary wave solutions whose orbits are confined 
inside the cells. The Picard Theorem about the existence and uniqueness of solutions of ordinary 
differential equations guarantees that there are no solutions of the analogous mechanical system 
crossing the edges of the cells, except through their vertices, where the theorem is not applicable, 
see (|3U|) . Therefore, the existence of singular kink orbits on the reticulum rules that any other kink 
orbit must be confined inside the cells in such a way that different orbits can meet only at the 
vertices of the cell. 



£[K BF + VF + F-VF-A ]M 



ri+rn 



n M*) a - 

n+m 

II Hxf-afr* 

n+m 

n n x f-air< 



3.2 Generic kink varieties 



Full knowledge of the kink manifold in this type of model is achieved by finding all the separatrix 
trajectories (finite mechanical action) of the Newton equations (|16|) equivalent to the static field 
equations. The most effective method of solving (|16p is the Hamilton-Jacobi procedure because 
these models are chosen in such a way that elliptic coordinates allow Hamilton-Jacobi separability. 
From the mechanical Lagrangian 



1,92 



l u 2 -v 2 (du\ 2 lvP_ 

2u 2 -n 2 \dx J + 2n 2 



, f(u)+g(v) 
dx J u 2 — v 2 



(33) 



the generalized momenta are determined: 



Pu 



dC 



u 



v du 



i 2 - n 2 dx 



Pv 



dC 

dm 



u 



n 2 



V 



2 dx 
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The mechanical Hamiltonian H = £n^i reads: 

H = -^— 5 {h u + h v ); K = \{u 2 -n 2 )p 2 u - f{u) , h v = \{tf -v 2 )p 2 v -g(v) 
u z — v z 2 2 

The Hamilton-Jacobi PDE equation 

ax \ou ov J 

corresponding to this Hamiltonian is separable. Denoting v! = ^ , v' = the separation ansatz 
i7 = <Jx(x) + J u (u) + Jv{v) provides solutions of the PDE (fM|) in terms of the solutions of the 
ODE's: 



/ is 2(F + Eu 2 + f(u)) „, . 2(-F - Ev 2 + g(v)) 

(35) 

because ([Ml) becomes: 

\{u 2 - n 2 ) {Xf - En 2 - /(«) = F = - v 2 ) {j'f - Ev 2 + g(v) . 

E is thus the mechanical energy; the value of the first invariant I\ = E. The relation of the 
separation constant F with the two first integrals is a little bit more involved. From 

h = Cnh= 2 l 2 {(u 2 - n 2 )h v - (n 2 - v 2 )h u } ; h u = F + Eu 2 , h v = -F - Ev 2 

one deduces that F is a linear combination of the two invariants: F = — Q 2 Ii — I 2 . 

The first equation in (|35p is integrated immediately to find J x = —Ex whereas the solution of 
the other two equations in ([35]) are the quadratures: 



f >\ /* , l2(F + Eu 2 + f(u)) „ . f , l2(-F-Ev 2 + g(v)) 
J u = sign («') y ^ 1 U 2_w ' ^ = Slgn M 7 * V ^ _ ,,2 ■ 

" " ( (36) 

Hamilton-Jacobi theory prescribes that "particle" orbits are determined by (j37[) (left) and "time 
schedules" ruled by ([3"7]) (right): 

|f =71 , % = 12 , (37) 

where ,y is the Hamilton principal function, the solution of (j34p and (|35p . and 71, 72, are integration 
constants. 

Finite action trajectories require that E = F = and give the solitary wave solutions, finite 
energy, in the field theoretical model. Therefore, the kink orbits, corresponding to finite action 
trajectories, and the kink profiles, given by the time schedules of these orbits, are the quadratures: 

sign(u') I ■= = d " = = - sign(v') [ = = ^ = = = v^i (38) 



y/(u 2 - n 2 )f(u) 'J ^(n 2 -v 2 ) g (v 

f u 2 du . . f v 2 dv f- s 

sign(n } I Tw^mm " Slgn(v ] J vw^rn = + 72) ■ (39) 

Thus, 72 obeys time translation and is associated to the first invariant. In the field theoretical 
context this parameter sets the center of the kink and/or the lump of energy. The meaning of 
71 is more subtle: orbits for different 71 are related to each other by non-linear transformations 
generated by the second invariant I 2 - 
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3.2.1 Linear stability analysis 



The qualitative structure of the moduli space of kink solutions (|38p - (|39p parametrized by the (71, 72) 
constants is explained by a linear (or local) stability analysis. We write the first-order ODE system 
(|30p in the form: 

^ = (-irF(u,v) , %L = (-l)f>G(u,v) . (40) 

Note that (—1)° = sign(u') and (— 1)^ = sign(u'). The fixed points in the trajectories are identified 
as the zeroes of the functions F and G: F(uq,Vq) = = G(uq, vq). Thus, the points in M. are fixed 
points and one can study the behavior of the trajectories near those points by linearization of the 
ODE system fl3D]): 



(41) 

It is easily checked that, depending on the choice of a and /3, the fixed points belonging to Ai fall 
into three categories: 1) stable nodes (the two eigenvalues of M^ a '^ are different and negative), 2) 
unstable nodes (distinct eigenvalues, both positive) 3) saddle points (one positive , one negative 
eigenvalue). The trajectories flow inwards and end in the fixed point (case 1), start in the fixed 
point and flow outwards (case 2), and flow first inwards, skip the fixed point, and go outwards (or 
viceversa, case 3). 

Moreover, taking the quotient of the two equations in (|40p . one obtains the differential equation 



du = (-i)°f( ¥ ) = (-ir Hv?-w)f{u) 

dv {-l)PG(u,v) (-1)^ {n 2 - v 2 )g(v) ' 1 ' 

showing that the flow is undefined, ^ |( UOil)0 ) = §, precisely at the fixed points. Thus, there are 
pencils of trajectories parametrized by 71 flowing in, flowing out, or skipping the fixed points. 

To distinguish which points of Ai are nodes or saddle points one must focus on the values of 
a and /3. If a = /3, the flow induced by ()42|) runs along monotonic increasing curves in the elliptic 
strip E b ecause ^jj is always positive. On applying the linear stability analysis to a type 1 cell, 
Figure 7(a), one finds that C and B are alternatively stable or unstable nodes if a = f3 = or 
a = (3 = —1 (red orbits in Figure 8(a)) whereas D and A are saddle points. The roles of C,B and 
A, D are exchanged (blue orbits) if a ^ /3. Given the existence of solutions on the borders of the 
cell, Picard's theorem confines inside the cell the pencils of kink orbits connecting the fixed points 
either on the vertices C and B or on D and A, depending on the relative values of a and /3, see 
Figure 8(a). These two families of kink orbits are parametrized by the integration constant 71 E M. 



3.2.2 Global stability and the boundary of the moduli space of kinks 

The stability of the solitary waves is inherited from the global -rather than local- stability of the 
kink orbits. By global stability we mean the stability of the trajectory as a whole, not merely the 
identification of the character of the fixed points along the trajectory. According to Morse Theory, 
see [jQ, [28], [21], [3H]: global stability (or instability) is essentially accounted for by the number of 
conjugate or focal points crossed by a pencil of trajectories between the starting and ending points. 
Conjugate points are the points where trajectories forming a congruence meet. Except in the Type 
I-2a and Type I-2b models, where the foci belong to Ai and are fixed points, the foci are not fixed 
points of ([30]) . The flow ([42]) at the foci is always undefined: ^ |( uo= i i „ 0== |-i) = §• In sum, in 
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Type I-la and Type I-lb models the foci are conjugate points to the stable and unstable nodes 
on the other vertices of the cell, and families of globally unstable trajectories cross these points. 
Therefore, the blue and red kink orbits in Figure 8(a) are stable, whereas the blue orbits in Figure 
8(b) and the blue and red orbits in Figure 8(c) are unstable. 




Figure 8: Families of kink orbits in each type of cell: a) K^ip/i) and K^ (71) in Type 1 cells, b) (71) 
and K^ A (ji) in Type 2 cells, and c) K% (jx) and K^iji) in Type 3 cells. 



Bearing this in mind, we shall describe the families of kink orbits inside the cells. In the cells 
where one or two of the foci are vertices that do not belong to M , kink orbits flowing according 
to (|42|) and connecting two elements of Ai are formed by gluing two families of solutions that 
correspond to two different choices of a relative to ft. The gluing must be performed in such a way 
that the kink orbits back in the M. 2 Cartesian plane are continuous. We now describe the different 
types of kink orbits arising in different types of cells. 

• In type 1 cells there are two families of kink orbits connecting opposite vertices in the 
cells, distinguished by sign((— l) a+/3 ) on (j4~2|) . We shall refer to these two types respec- 
tively as i^ D (7i) when sign((-l) Q+ 0) = -1, and i^ 5 (71 ) , if sign((-l) a+ ^) = +1. When 
sign((— 1) Q+ ^) = —1 A is a stable node, and D is an unstable node, or viceversa. In this case 
B and C are saddle points. If sign((— 1) 0+/3 ) = +1 the pair of points A-D exchanges its role 
with the B-C couple, see Figure 8(a). The parameter 71 labels one member in each of the 
two families of kink orbits. The subscript suggests that these kinks are formed by two basic 
kinks/lumps, which can be confirmed by taking the asymptotic values: 71 — > ±00. 

To prove this last statement, we consider a solution characterized by the value 71. Therefore 
for generalized MSTB models (|38p the orbit of this solution complies with: 



/^a f ^ /_ 1 n / 8 f ^ _ ^2 

J (u 2 -SlV»nf=i l« 2 "^h J {^l 2 -v 2 )\v\^\[^ =1 \v 2 -a}\^~ V " 

(43) 

Linearization of (|43p respectively near the point C reads: 
where we have defined: 

1 

C u (ai, . . . , cr n + m ) 



C„(cri, . . . , a 



1 



5 O n + m ) — 



Likewise, linearization of (|43p respectively near the point B reads: 

V271 = (-1)° / du B ^---' an+m) -{-if I dv B ^U---^n+m) (4 ) 
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where now we have 

B u {a,,...,a n+m ) = — - n fc_i w+i _ <T 2 )ai(ff 2 +i _ a 2 )Qfc ng^| _ al +1 )^(2a k+1 r^ 

B v {ar,...,a n+m ) = — _ fjgT^ _ _ ^ flS^K - ^f^R^j^ ' 

In (|4"4"j) and (|4"5|) . C u > 0, C v > 0, 5 U > 0,B V > and the pertinent integrals are simple 
although we must distinguish between the cases a r = 1 or a r > 1 with r = j,j + 1, fc, k + 1. 
Whichever the case, the limit 71 — > 00 in (|44p and (|45p provides us, for instance in the case 
a = /3 = 1, with the orbit v = a; in (|44p . and u = Ofc+i in ()45p . while the limit 71 — )• —00 
yields the orbit v = (Tj+i in (|45|) and u = in (|44|) . This means that the combination 

R BD + R DC Qf 

singular kinks arises as 71 — ^ —00 in the moduli space of K^ G kinks while 
the combination K.f A + K AC lives at the 71 = 00 boundary of the K§ kink moduli space. 
The same analysis can be tediously repeated in the rest of the cases. In sum, we find that: 

The 71 parameter, giving the orbit in the mechanical problem, measures the distance between 
the centers in the combination of basic kinks. Finite values of the parameter characterize the 
nondinear combinations of single kinks as the finite separation between lump centers. The 
kink families K^^i) and K^iji) are formed by different compositions of two kinks at 
different distances that become the basic or singular kinks when the inter-center distance is 
infinite. 

• In type 2 cells there are also two families of kink orbits but one of the families is of a new type. 
If sign((— = +1, a K^iji) family arises with identical properties to the K2 (71) 
family in type 1 cells. B and C are respectively stable and unstable nodes, or viceversa, and 
A is a saddle point, see Figure 8(b). 

The novelty comes from the fact that in this type of cell, D is a focus, F, that does not 
belong to M because it is not a fixed point of the trajectory flow. F not being a fixed point, 
neither the K AF {"ii) nor the K^^i) kink orbits -expected when sign((— 1) Q+ ^) = —1 and 
A is either a stable or unstable node- can exist. In Figure 8(b) the apparent homoclinic 
kink trajectories K AA {^f\) are plotted in blue. These trajectories correspond to one K AF {^fi) 
and one K F (71) kink orbit glued at the focus. One might think that such orbits fail to be 
derivable at the focus and consequently that their action is infinite, yielding field theoretical 
solutions of infinite energy. To examine this issue closely, it is convenient to go back to the 
Cartesian plane. 

Recall that the inverse map from E 2 to M? induced by the change to elliptic coordinates is one- 
to-two, see (f20|) ): E 2 is mapped into the upper half-plane <p2 > by £ + * and also into the lower 
half-plane 02 < by <!;~*. Therefore, the inverse image of the point A = (<r n _i, er n +i) is the 

two-element set: A ± = ( ^-i^+i ; ± V^E^^E^ Ay The ^(71) can be chosen to be 
mapped in the upper half-plane of the Cartesian plane M 2 or viceversa. These two choices give 
rise to a family of continuous and derivable kink trajectories connecting the stable node A + 
with the unstable node A~ (or viceversa when A~ is stable and A + is unstable) H. Therefore, 
every kink trajectory of this type is of finite energy, heteroclinic despite of appearing to be 
homoclinic in E 2 , and crosses from the upper half-plane to the lower half plane or viceversa 

5 If the two K^ F , K FA families of the connected sum Kt A {-y\) = K%* '(71) \j K FA (^\) are mapped into the same 
half-plane of R 2 , the kink trajectories fail to be derivable at the focus and develop infinite action. 
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through a single point: the focus F. The last point confirms that F is a conjugate point to 
both A + and A~ and, as anticipated, these trajectories are unstable according to the Morse 
index theorem. A laborious analysis, but conceptually identical to that shown for type 1 cells, 
shows that at the | -yi | — > oo limit the following kink combination arises: 

lim Kt A { ll ) = Kt C + K^ FB + K^ A 

|7i|->oc 

C F \l F B 

where by K 2 we denote the two-step singular kink orbit living respectively at the edges 

u = 0, and v = in two different periods of time, the second period synchronized with the 
previous one at the focus F = (O, fi). Note that 

r^CF\/FB ,. tsCB( \ 

K 2 = hm K 2 (71) 

|7i|->oo 

and the two last limits above account for the subscript 4 in the nomenclature. In R 2 , the 
j^CF\J fb jjj^ qj.]-,^ \[ ves on the 02 = axis. 

The translation from this analysis in the analogous mechanical system to the field theory is as 
follows: 1) To this one-parametric family of finite-action heteroclinic kink orbits corresponds 
a one-parametric family of finite-energy topological kinks. 2) Because the kink orbits are 
unstable, the topological kinks are also unstable. 3) These kinks are composed of four basic 
kinks. 

In type 3 cells only two of the four vertices, A = (— f2, 02), B = (SI, 02), are elements of M 
whereas the other two vertices are the foci, F± = (±0,0), of the elliptic coordinate system, 
that do not belong to M. If sign((-l) a+/3 ) = +1, A is a stable or unstable node whereas B 
is a saddle point. When sign((— l) a+l3 ) = —1 the roles of A and B are exchanged. Therefore, 
except for singular kink orbits all the generic kink trajectories cross one of the two foci, which 
are not fixed points but display indefinite flow. 

For these reasons, the features of the kink trajectories are better understood in the Cartesian 
plane. In R 2 the vertices of the cell become: A = (—02; 0), B = (02,0), and F± = (±0,0). 
Kink trajectories in the K^iji) blue family start from A, run along the upper or the lower 
half-plane, meet at the F + focus, where they pass to the other half-plane, and come back to 
B. This tour is performed either clock- or counter-clockwise. The red kink trajectories in 
Figure 8(c) are identical to the blue kink trajectories of the same figure, but B is the starting 
and ending point, and the bridge between half-planes is the F_ focus. The kink orbits in type 
3 cells are thus unstable homoclinic trajectories: each trajectory starts from either A or B 
and ends at either A or B. The foci F + and F- are respectively conjugate points to A and 
B where all the members of the family meet. Again, it is hard but merely routine work to 
prove that the |7i| = 00 limit is the following combination of singular kinks: 

lim K£ A ( ll ) = K? B + Kf-" JF - F+ " JF+B 

|7l|-K» 

lim K* B { ll ) = K* A + K BF +V F + F -" JF - A 
|7i|-K» 

j^AF-\J F-F + \J f+b ^gj^gg three-step singular kink orbit living respectively at the edges 

v = — fi, u = O, and v = O in three different periods of time, the second period synchronized 

with the previous and later periods at the foci. Until now we have met only kinks for which 

the number of the steps coincides with the number of lumps. For instance, the kink orbit 
CF \/ F B 

K 2 runs over two edges, and even though the energy density of these stable kinks is 
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formed by a single lump, it is exactly the superposition of one and one Ky kinks 

A.F \I F F \l F B 

with coinciding centers. K 1 - + v + kinks differ in two ways: 1) Combinations of 
three K^ B do not form these kinks. 2) K^ F ~ V F - F +\l F + B ki n k s are unstable because upon 
crossing one of the two foci they pass a conjugate point of the congruence of kink orbits they 
belong to. In sum, there are no reasons for this configuration to be considered as a composite 
kink and, moreover, Kf F ~ * F ~ F+ V F+B kinks are unstable. 

The translation from this analysis in the analogous mechanical system to the field theory is as 
follows: 1) To this one-parametric family of finite action homoclinic kink orbits corresponds a 
one-parametric family of finite energy non-topological kinkg^|. 2) Because the kink orbits are 
unstable the non-topological kinks are unstable. 3) These kinks are composed of one basic 
kink and one unstable singular kink, see 

The previous results can be summarized as follows: 

Proposition 3. There are families of kink orbits confined in each cell in the elliptic plane. The 
kink trajectories connect elements of Ai at opposite cell vertices. We shall denote these generic kink 
orbits by labelling the vertices connected by them: K^iji), K BC (^/i). There are two special cases: 
either one or two cell vertices are not zeroes of V(</>i, fa). In type I-la and type I-lb generalized 
MSTB models type 2 cells arise: in these cells there is a family of topological kink orbits connecting 
the elements A + = a n+ \) and A~ = £~(er n _i , o- n +i) of M. in M 2 ; every orbit in this family 

passes through one of the two foci. In n = 1 type I-la generalized MSTB models there is one cell of 
type 3; in this cell there is a family of non-topological kink orbits joining the points A = (— 0,02) 
or B = (f2,o"2) with themselves. Each non-topological kink orbit crosses one of the two F± foci. 



3.3 Sum rules of kink masses 

The action of the mechanical trajectory is the energy of the static field profile. Therefore, 



E\ 



. 2(F + Eu 2 + f(u)) 
du sign(n );/ u 2 _ n 2 + 



+ 



*4> 



, 2(-F - Ev 2 + g(v)) 



ft 2 



which is the Hamilton principal function -the action- of the mechanical trajectories, is also the en- 
ergy of the static field solutions. Solitary wave solutions have finite energy E[4>\, fa] and correspond 
to trajectories with zero mechanical energy and zero separation constant: E = F = 0. 

Proposition 4. The energy of the kinks in generalized MSTB models is the mechanical action 
of the kink orbits: 
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>l,92 



du 



2/(«) 



u 2 -O 2 



+ 



dv 

TTi, (£*</>) 



2gM 



n 2 



(46) 



where ir u (^*(j)) and tt v (^*4>) are respectively the projections of the point of the kink orbit into 
the axes u and v. 

Therefore, all the members of the kink family are iso-energetic. Formula (|46|) allows us to state 
the following: 



6 Amazingly, the non-topological kinks of the MSTB model were discovered at the very beginning of the search for 
solitary waves in field theories with two real fields. Of all possible kinks that we have just described these are by far 
the most bizarre. 
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1. The kink energy does not depend on the fine details of the solitary waves: kinks with very 
different profiles share the same energy. The kink energy is determined only by the projections 
of the kink trajectories into the axes. 

2. The important implication of this is that the energy of any member of a family of composite 
kinks is uniquely determined as the sum of the energies of the basic or singular kinks. The 
total energy of the singular kinks is: 



E[K AB ] = E[K[ 
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/ dv 



2fl(«) 



n 2 



E[K AC ] = E[K( 



BDi 



CT fc+i 

du 



2/(«) 



n 2 



(47) 



A little more care is needed when the foci are crossed by a congruence of kink trajectories. 
Here, line integrals must be computed piecewise: from the starting point to the focus and 
from the focus to the endpoint. For instance, the energy of K^ F ^ FB kinks in Type I-la and 
I- lb generalized MSTB models is: 



E[K. 



CF V FBt, 
2 J 
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dv 



2g(v) 

n 2 -v 2 
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(48) 



Simili modo, the energy of unstable singular kinks Kf F ^ ' F F+ ^ F+A in type 3 cells of n = 1 
Type I-la generalized MSTB models is also computed piecewise: 



E[K AF-\,F.F + VF + 
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2g(v) 
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(49) 



Kinks or combinations of kinks that have the same projections on the elliptic axes have the 
same energy. This is the crux of the matter of the kink mass sum rules. We now write all the 
kink mass sum rules observed in the generalized MSTB models. These sum rules connect the 
energy of the kinks solutions confined in a cell, and hence we must distinguish three possible 
cases: 



Kinks confined in cells of Type 1: 
E[K? B ] = E[K^ D ] 



E[K? C ) 



E[K{ 



BD-i 



E[K£ D ( 7l )] = E[Ki c ( 7l )} = E[K? B ] + E[K? C ] 



Kinks confined in cells of Type 2: 

E[K^ FyFB ] = E[Ki B ] + E[K? C ] 



E[Kt A {li)] 



. 4 y- nji -2E[K^]+2E[K( j 
Kinks confined in cells of Type 3: 

'AA(„,.x\ _ vlisABi , v]T sAF-\J F-F+M F+B, = E BB 



E[K* A { ll )] = E[K? B ] + E[K? 
E BB ( 71 ) = E[K BA ] + E[K B F+ VF+F " VF - 



(71) 



These kink mass sum rules are due to the Hamilton-Jacobi separability in elliptic coordi- 
nates of the mechanical systems ruling finite energy static solutions in the generalized MSTB 
models. 



3.4 Kink orbits as gradient flow lines 

It is interesting to address the results of the previous sub-Section from the point of view sketched 
in the Introduction. The Bogomolny'i-Prasad-Sommerfield understanding of several types of topo- 
logical defects in field theory, such as kinks, special classes of vortices and magnetic monopoles, 
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and instantons, can be applied to some solitary waves of generalized MSTB models. We recall that 
the BPS approach is possible if the potential energy density is equal to half the norm of the square 
of the norm of a function W(4>i, <fe), see ([2]): 



u(<k,<h) 



1 fdW dW dW dW 



+ 



2 \ucpi 

Writing the static energy in the form ( 
E$\ = \ 

one sees that solutions of the first-order ODE system 
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(50) 



(51) 



dW _ 



(52) 



are absolute minima of the energy if W behaves well enough. Stable kink orbits are thus the 
flow lines induced by the gradient of W0, usually referred to as the superpotential in physicists' 
literature, because of the central role that this function plays in supersymmetric models. A closer 
look to (|50p reveals that W is no more than a solution of the time-independent Hamilton-Jacobi 
equation of the analogous mechanical system with zero mechanical energy: i.e., the superpotential 
is the E = Hamilton characteristic function in Hamiltonian dynamical systems. 

For mechanical potential energies of Liouville Type I systems, (|50p in elliptic coordinates reads: 



l(u 2 -v 2 )\ y '\ 
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(53) 



The separation ansatz W(u,v) = W\{u) + W2(v) reduces the PDE 
order differential equations 



to two independent first- 



dWi 



2f(u) 

du v ' V u 2 — ft 2 
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dW 2 
dv 
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2g(v) 
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which can be integrated by quadratures: 



w {a '^{ u ,v) = (-iy 



du 
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(54) 



Proposition 5: Generalized MSTB models are built from the complete integral - J5^| )- of the 
time- independent Hamilton-Jacobi equation with zero mechanical energy for the analogous mechan- 
ical system. In modern language, we characterize the generalized MSTB model as those admitting 
four different superpotentials according to the freedom of choice of signs in (5$ . 

More precisely, 

Proposition 6: In E 2 the four superpotentials of generalized MSTB models have the form: 



C n W = W 1 (u) + W 2 (v) 



d 2 (C Q W) 
du dv 







(55) 



7 The minima of U, the elements of M, are either minima, maxima, or saddle points of W. Accordingly, the 
critical points of W are either stable nodes, unstable nodes, or saddle points of the gradient flow. 
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In M. 2 the PDE A55\) becomes the cumbersome PDE satisfied by generalized MSTB superpotentials: 
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(56) 



The concept of superpotential allows us to write the invariants (|16|) of the mechanical system in 
the surprisingly simpler form: 
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(57) 
(58) 



explicitly showing that I\ = = Ii for the solutions of ([52]) , 

Obviously, because four superpotentials are available there are four different ways of saturating 
the BPS bound [7]. For instance, in E 2 these equations are: 



d_u = ^/2{v?-W)f(u) 



dx y ' u 2 — v 2 



dv = , l)0 ^/2{Sl 2 -v 2 )g(v 



dx 



u 2 — v 2 



(59) 



It is clear that for kink orbit solutions of (|54p the kink energy 



dW 



K 



the integral of the exterior derivative of the superpotential along the kink trajectories, is precisely 
given by formulas (|46p . ()47p . (|48[) . and (|49|) for the different types of kinks. If the kink orbit does 
not cross points where W is not differentiable Stoke's theorem rules that 



E\ 



and the energy of stable kinks is a topological charge. The situation is more subtle for kink orbits 
that cross non-fixed conjugate points. Stoke's theorem must be applied piece-wise, and we find: 



E[(f> K ] = W((t>i) - W(TOei)) + Wfrftei) - W(<j> 



respectively for unstable kink orbits in Type 2 and Type 3 cells. The important point is that these 
kink orbits are not solutions of (|59p for the same choice of a and /3 over the whole real line BL 
The unstable trajectories solve ([54]) for a given f3 relative to a between x = — oo and x = xo, 
the point where the focus is reached. Between xo and x = +oo, the first-order equations satisfied 
correspond to the other relation between f3 and a. The unstable kink orbits, however, satisfies the 
static second-order equations (|16p but the unstable kinks do not saturate the BPS bound. 

We show explicitly that the existence of four superpotentials is connected to the Hamilton- Jacobi 
separability of the analogous mechanical system to generalized MSTB models in elliptic coordinates. 
Formulas (|26p and (|27p for the two invariants in elliptic coordinates due to HJ separability provide 
the following identities for the square of u and v derivatives: 
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The separatrix trajectories between bound and unbound motion arise for I\ = 0, I2 = 0, values for 
which the above equations reduce to: 



du\ 2 _ 2f(u)(u 2 -Q 2 ) fdv\ 2 _ 2g(v)(n 2 -v 2 ) 

dx) (u 2 -v 2 ) 2 ' \dx) (u 2 -v 2 ) 2 ' 1 ' 

equivalent to (|59p with the four sign combinations. Complete Arnold-Liouville integrability is 
accompanied by HJ separability, allowing us to obtain all the separatrix trajectories - all the kink 
profiles- from the complete solution of the HJ equation: the four superpotentials. The partial 
solvability of the HJ equation would permit only a partial identification of the kink manifold. The 
complete solution of the HJ equation for arbitrary values (non-zero) of the invariants, 



w^> = nr / iu^Eh ± SgEgS ± /(-I J du . /*»+/.<!«" -»>)+<*.)) 



n 2 -v 2 



also reduces to algebraic equations between quadratures all the other trajectories - either periodic 
or unbounded- of the analogous mechanical system. Back in field theory, the periodic orbits are 
interpreted as (unstable) kinks on a circle (the space-time being ExS 1 with hyperbolic metric) 
whereas unbound mechanical motion corresponds to static field solutions of infinite energy. 

To finish this Section and a thorough study of generalized MSTB models we now address this 
issue in M. 2 , i.e., we analyze the I\ = = 1% condition in Cartesian coordinates. The surprise is 
that not only it is satisfied by solutions of the first-order ODE system ([52]) . but that solutions of 
the new first-order ODE system 



dx ^—^f+mK+w+ A2 
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(61) 

2J 



also comply with I\ = = /2- This implies the existence of a second superpotential W. The 
vector field from the right hand sides of (|6ip : 

A(<t>U<l>2) = =F 



- n - T 



V[(^-n) 2 + ^][(^i + n) 2 + ^] A /[(0 1 -r2) 2 + 2 ][(0i + o) 2 + ^ 



despite its horrible aspect, is curl-less because equation ([56]) is valid in Type I Liouville models. 
Therefore, Green's theorem ensures that there exists a function W such that (|6ip is equivalent to 
the first-order ODE system: ^ = |J, ^ = §£• Four different superpotentials, ±W, ±W, are 
also available in Cartesian coordinates. 



4 Eighth-order generalized MSTB models 

The general analysis of the previous Section affords us a practical method of obtaining the kink 
variety of generalized MSTB models by means of a few computations. The references [16-28] are 
devoted to the original MSTB model, in which the potential energy density is a quartic polynomial 
in the two scalar fields. In the paradigm of generalized MSTB models, the functions / and g 
are respectively: f(u) = \{u 2 — l) 2 (u 2 — fi 2 ) and g{v) = \{v 2 — 1) 2 (0 2 — v 2 ). Therefore, M g = 
{er_i = — VL;o\ = 0} and A4f = = 0; 02 = 1} such that m = 1 and n = 1, a type I-la 
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generalized MSTB model that involves only the type 3 cell Cl 2 1 in the reticulum TZ. There are 
only two singular kinks: K^ B and Kf F ~ ^ F ~ F+ V F + B (respectively denoted as TK2 and TK1 kinks 
in the literature). The remaining kinks form a family of non-top ological (unstable) kinks: K^ A 
(usually referred to as NTK2 kinks). In other papers, [29-31], a generalized MSTB model where the 
potential energy density is a degree-six polynomial in the fields is discussed. The / and g functions 
in this model are: f(u) = \u 2 (u 2 — l) 2 (u 2 — a 2 ) and g(v) = \v 2 (v 2 — l) 2 (d 2 — v 2 ). A type I-lb 
generalized MSTB model with m = 1 and n = 1 arises. The set of zeroes of g and / are respectively 



Mg = R 



-Q; (Tq = 0; o"i = $7} and M. 



j = {cTi=f2;cr2 = l} whereas the set of minima of the 
(±1,0);5± = (0,±VT^2);0 = (0,0)}. 



potential in R is, using the notation of [36J, M. = {A± 
The reticulum in R 2 is formed by two type 2 cells joined by the edge on the ordinate axis. There are 
three singular kinks, denoted as K^ , K BO and K± B , and two families of topological composite 
kinks, one of them stable, K^i'Ji), and the other one, K BB (~{i), unstable. 

To close this article we shall study a two-parametric family of generalized MSTB models with 
an eight-order polynomial as potential energy density: 
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where a = ^f(2Q 4 - 2a 2 Q 2 - 2t 2 Q 2 + 4fi 2 - a 4 



r 4 ) and Q 2 = a 2 r 2 . 





Figure 9: Potential density energy in the model (t < a < 1) and bifurcation of the elements of M.. 

The parameters a and r are the coupling constants of these generalized MSTB models and the foci 
of the elliptic coordinate system are determined from their product. The four superpotentials are: 



w{<h.,<h 



:^|^i±^) 2 + 0| [3(01 + 4>2=f nfa + n 2 )^ 2 + <j> 2 2 ± + n 2 y 

-5(a 2 + r 2 )(0 2 + 2 =F O0i + n 2 ) =f 3ft0i(0 2 + (j>l + CI 2 ) + 15ft 2 ] 



■ (63) 



The secret to finding a complete solution to the zero mechanical energy Hamilton-Jacobi equation 
of the analogous mechanical system is that the mechanical system is of Stackel type. There are 
two invariants in involution -the mechanical energy, /]_, and a generalized momentum, 1^- that are 
quadratic in the momenta: 



h = 
h = 



dx 



d.r 



The mechanical system is amenable to the separation of variables. Here, /3 = ft 2 (ft 2 — a s 
I) " lO* + t% 5 = ft 2 - a 2 - t 2 and « = 2^±d. 
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In the count and localization of the elements of Ai a complicated process of bifurcation arises, 
depending on the values of the parameters r and a. To analyze the most interesting cases, we shall 
consider a fixed value of a £ (0, 1) and let allow r to vary along the real line. If both parameters are 
greater than or equal to one the model presents only isolated singular kinks. There are three cases 
- 1) t< 1; 2) t = 1, and 3) r > 1- where degenerate kink families arise in the model. In Figure 9b 
we show the elements of At when r is varied from to 3. If r E (0, 1), At is formed by four elements 
living at the intersection points of the two enveloping ellipses of the separatrix trajectories with the 
0i-axis. For r > 1 there are still four elements of Ai in the </>i-axis but also another four elements 
of Ai at the intersection of the enveloping ellipse and hyperbola of the separatrix trajectories arise, 
r = 1 marks a pitchfork bifurcation of the elements in the interior ellipse when this curve becomes 
a hyperbola for r > 1 . All these statements become clear looking at the potential energy density 
in elliptic variables: 

CrUfaih) = [(« 2 ~ - 2 ) V - r 2 ) V - + {v 2 - V - r 2 )\tf - v 2 )] (66) 

Comparison between ([BU]) and ([25]) shows that we are dealing with a family of generalized MSTB 
models. 

4.1 The structure of the set Ai 

We shall now describe the distribution of the zeroes of the potential energy density both in E 2 and 
M 2 successively for r < 1, r > 1, and r = 1. 

• Regime I: < r 2 < a 2 < 1. The model belongs to type I-la generalized MSTB models. 
Because Q, < t < a we have: n = 1 , m = 2. The set of roots of g and / are respectively 
At g = {(T-i = — 17; o\ = £1} and Aif = {ai = $7; a 2 = r; 03 = a}. Therefore, the set of zeroes 
of U(c/>i, $2) has four elements, both in E 2 and M 2 , in this range of the parameters: 

Mi = {A ± = (a,±n)-B± = (T,±n)} , Mi = {A± = (±a,0);B± = (±T,0)} . 

The model exhibits symmetry with respect to the Klein four-group Z 2 generated by the 
reflections of the fields: <pi — > —<j>\ and $2 —> —<p2- In the quantum version of the model, the 
symmetry with respect to the Z 2 sub-group generated by <j>i — > —<pi, however, is broken by 
the choice of any A or B ground state. The set Aii is formed by the A- and S-orbits of this 
sub-group, such that the moduli space of vacua has two elements: JVlifL 2 , = {A, B}. 

In Figure 10 the reticulum corresponding to this regime is shown. It is formed by two cells: 
Cl 2 n , a type 3 cell delimited by the edges B + B-, B + F + , F + F^., and C 23 n , a type 1 

cell whose vertices are the points A + , A—, B + , and On the Cartesian plane, M 2 , the cell 
is the region confined between the two ellipses 

ti + § = r 2 (left) , 4>l + f 2 =° 2 ( ri ght) , (67) 

whereas is the region bounded by the ellipse (|67p (right). To simplify this and subsequent 
formulas we use the following conventions: a = yl — a 2 , f = \/\\ — r 2 |. 

• Regime II: < a 2 < 1 < r 2 . Because o < fl < r, the model also belongs to Type I la 
generalized MSTB models in this regime, but now: n = 2, m = 1. The set of roots of g and 
/ are respectively: M. g = {o\_2 = — SI; cr_i = — o;a\ = o\G2 = -M-f = {02 = J7;o"3 = r}. 
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The cardinal of the set of zeroes of the potential energy in the elliptic strip E 2 differs from the 
same number in the Cartesian plane M. 2 . Mn has six elements whereas Mu contains eight: 

Mn = {A± = (r, ±Q) ; B± = (Q, ±a) ; C| = (r, ±o) } 

Ma = {A± = (±t,0); B± = (±a,0); Cf = (1,±0); C± = (-1,±0)} , fi = a|f | . 

The key point in the bifurcation process is that the ellipse (|67p (right) arising in regime I 
becomes a hyperbola in regime II. As in regime I, there are two orbits -one A, and one B- 
in .Mn of the Z2 sub-group generated by <j)\ — > —<j>\. The four C points, however, form a 
single orbit of the full Z| field reflection group because the C + and C~ orbits are connected 
by the other reflection <p2 — > —<p2, which would also be spontaneously broken in the process 
of quantization if one of the ground states, C, is chosen. 

Figure 11 shows the reticulum in E 2 in this regime. It is formed by three cells: C 23 n , C 2 2_i, 
and C 2 |. The images of C^-i an d C 2 | in R 2 are connected by the field reflection <pi — > —cfri. 
C±2±i are type 2 cells with vertices at the points A±, B±, C± = C±, and F±. The vertices 
of the type 1 cell C 2 ^ are the points B-, B + , C_, and C¥. 




Figure 11: Reticulum in E 2 and elements of M. in the regime II. 

• Regime III: < a 2 < 1 = r 2 . This is the border between the two previous regimes. The 
potential energy density ([66]) becomes 

CU= — L-j [{u 2 - a 2 f{u 2 - l) 2 - (v 2 - a 2 f(v 2 - l) 2 ] . (68) 

Because a = £1 < 1, the model belongs to type I-2a generalized MSTB models, with n = 1, 
m = 1. There are four elements in both M and M: 

M m = {A ± = {l,±n); B± = (n,±n)} , Mm = {A ± = (±1,0); B± = (±fi,0)} 

The phase transition occurs when the points B cross the foci, which are zeroes of the potential 
U. In Figure 12, the reticulum in E 2 for this regime is shown. There is only the type 1 cell 
C-ll) with vertices at the points A + , A_, B + , and B- 
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4.2 Kink Variety in Regime I 
Singular kinks 

Let us try the orbit 02 = in the second-order static field equations ([To]) . A simple integration plus 
some technical manipulations provide several basic kink solutions encompassed in the formula: 



0" + 01 
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T + (f>l 



,±2V2(<t 2 -t 2 )x 



(69) 



1(a). K BF V FF V FB : if |0 1 (x)| < r for all finite x and lim^-oo 0i(x) = ±r, \mi x ^. +00 0i(a;) = 
=Fr, the solution in (j69j) describes kinks K 1 that join the points B + and -B_, crossing 

the foci F± such that 0i = 0, -where the field changes sign- occurs for x = 0. These kinks are 
analogous to the TK1 kinks of the MSTB model. Their energy is: 



e[k bf\i FFy FB] 



\w(-t,o) - w(-n,o)\ + |w(-n,o) - w(n,o)| + \w(n,o) - w(r,i 
4V2 



15 



r[5ft 2 (l + 3a - a A ) - r 4 + fiV(3<7 d - 5a)] 



Because all of them cross a conjugate point (a focus), TK1 kinks are unstable. 

1(b). K± ±B± : if t < |0i (2;) I < a for all finite x and lim. c _ 5 ._ 00 0i(x) = ±r, lim x _^ +00 <f>\{x) = ±a 
the solution in (|69p describes K^ B kinks. The if^-kink orbits join the points A± with B± along 
the 0i-axis , see Figure 10. Their energy is: 



E[Kf ±B± ] = \W{±a,0) - W{±t, 



2^2 
"iT 



a -r) 3 (a 2 + 30 + r 2 



These kinks are stable because their orbits are not crossed by any other kink orbits. 

J^. _____ 

2. K l T : plugging the elliptic orbit (|67p (left) into the equations (fT6j) . solving in favor of 0i, 
and integrating the quadrature, we obtain: 



T + 01 



1 + 



(70) 



The K BB solitary waves characterized by this equation have kink orbits joining the B± with the 
Bzp points through the upper and lower half-ellipses (|67|) (left), see Figure 10. The second field 
component is given in terms of the first as: 02 (%) = ±<T\/ t 2 — 2 (x). The energy of the K BB kinks 
is: 



E[K[ 



B±B^ 
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3. K l ± T : by trying the elliptic orbit (j67[) (right) in the equations (|16p. solving in favor of 0i, 
and integrating the quadrature, we find: 

<j + 4>i a 

a-<f>\ 

The K AA solitary waves determined by this equation have kink orbits joining the A± with the 
A T points through the upper and lower half-ellipses ([67|) (right), see Figure 10. The second field 
component is given in terms of the first as: 4>2{%) = zt-fy 7 a 2 — 4>±(x). The energy of the K AA kinks 
is: 

E[K* ±At ] = \W(±a, 0) - W(T<r, °)l = ^HT n ^ 15 ~ ^ ~ ^ + ^ 

Amazingly, these kinks have the same energy as the K BB kinks despite the complicated expression 
of the super-potential W (|63|) . Two observations help to understand this point: (1) The expression 
(|63p for W is invariant under the exchange df>T. (2) The kink orbits K AA and K BB have exactly 
the same length in E 2 , see Figure 10 (right). 

Finally, the K AA and K BB kink orbits are not crossed by other trajectories. Thus, these kinks 
are stable. 
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Generic kinks 



By applying the Hamilton-Jacobi procedure to find the kink orbits and profiles respectively as the 
quadratures of (|38p and Q39|). we identify the kinks confined inside the cells. In this case, integration 
of the kink orbit equation (|3"8]) corresponds to 
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U + T 



U — T 



2 (<t 2 —t' 2 ) ( (J — U \ t 2 (ct 2 -t 2 ) 
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Q, — v \ ct 2 t 2 / r + v \ ct 2 (<t 2 -t 2 ) / a — v\ t 2 (<t 2 -t 2 ) 



Sign(w') 
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T — V 



a + v 



71 



(72) 



whereas the spatial distribution of these solitary waves, the result of integrating (|39p . is given by: 



'u — 17 \ S 2 T 2 



u + Q 



u + r 



t — u 



a 2 (a 2 — T 2 ) I a — U \ t 2 (<t 2 —t 2 ) 



a + u 



Q, — v \ * 2 t 2 /r + u\ <T 2 (<r 2 -r 2 ) f a — v\ t 2 (<j 2 -t 2 ) 



T — V 



a + v 



Sign(«') 
> X 

x Sign(v') 



(73) 



4. iv:f FVFiJ (7i): In the Ci 2 n cell — r < fa < r, -tct < cp 2 < ra in M 2 , or, fi < n < t, -fl < 
w < £7 in E 2 - there are no trajectories determined by (|72[) and (j73[) with either Sign u' = Sign u' or 
Sign u' 7^ Sign u' that connect nodal B points. Nevertheless, solutions of the second-order static 
field equations are obtained by continuously gluing one trajectory from a B point to the farther 
focus with another trajectory coming back from that focus to the B point. Necessarily, one piece of 
the joint two-step trajectory is determined by (|72|) and (|73|) with a combination of signs, whereas 
the second piece obeys the same equations with the other sign combination. 

These solitary wave solutions form a family of non-top ological kinks. There are two real inte- 
gration constants: 71 E (—00,00) and 72 E (—00,00). The second one sets the kink center, but 
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71 distinguishes between the members of the kink family. All these kink orbits connect the points 
B± with themselves, passing through one of the foci, see Figure 13. We shall refer to the members 
of this family as K BB (^i) kinks. The behavior of these solutions is completely analogous to that 
of the NTK kinks of the MSTB model: all the kink orbits meet at a single focus and are thus 
unstable. 




Figure 13: if 2 BS (7i) Orbits: in 
of cell are plotted in Figure 1. 



(left) in E 2 (right). The kink profiles and energy densities in this type 



Fi F \/ F B 

All the non-topological kinks are degenerated in energy, E[K 2 (71)] = | W(— r, 0) — W(—£l, 0)\-\ 

\w(-n, 0) - w(n, o)\+\w(n, o) - w(t, q)\+\w(±t, o) - w(tt, o)| = e\k? f yFF yFB ]+E[K B } = 

t 3 (10<t 2 + 45cr 3 — 15a 5 — 2t 2 — 15<rSl 2 + 9<7 3 S7 2 ), and their energy is equal to the sum of the ener- 
gies of the topological kinks living at the boundary of the cell. In keeping with this kink mass sum 
rule, one can check that these topological kinks arise at the 71 — > ±00 limit of the non-topological 
kink moduli space (j72|) : 

lim ^( 7 i) = Kf FVFFVFB + ^f s . 

71— >±oo 



5. K^ ±Bt (7i): In the cell C 23 n , the domain in M? bounded by the ellipses (f67|) fleft) and 

([67D(right), - 4>\ + || < r 2 , (fl + || < a 2 in M 2 , or, r < u < a , -0 < v < n in E 2 - there are 
curves characterized by (j72[) and (j73[) with the same sign combination joining the nodal A and B 
points. These solitary wave solutions form another one-parametric kink family. The kink orbits 
connect the points A± with B T or viceversa -depending on the sign combination- and are thus 
topological kinks, see Figure 14. The orbits are not crossed by any other trajectory and these kinks 
are stable. Their energy is: E[K^ B {-n)] = \W({±a,0) - W(^a,Q)\ + \W((±a,0) - W(±r,0)| = 
E[Kf A ] + E[K AB ] = |W((±r,0) - W{^t,0)\ + \W({±a,0) - W(±t,0)\ = E[K BB ) + E[Kf B ] = 
^^(cr 5 — bcrO, 2 + 5f2 2 r + 150 3 — r 5 — 5$7 3 r 2 + 3S1 5 ); i.e., equal to the sum of the energies of two 
of the topological kinks living in the cell boundary. Again, it is easy to check that these singular 




Figure 14: K^iji) kinks: Profiles and energy density (left) Orbits in IR 2 (center) in E 2 (right). 
kink combinations arise at the boundary of the K 2 B kink moduli space: 



lim K£ B { 1X ) = { 
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In sum, there are six different types of kink solutions: four single kinks, the K BF * FF V FB ; 
the K AB , the K AA and the Kf B kinks, and two kink families, the Iff- 5 (71) and K£ B {n) kink 
moduli spaces. The boundary of the K BB {^\) family is a non-linear combination of the single kinks 
j^bf\/ ff\/ fb x BB whereas the boundary of the K AB (ji) family is a non- linear combination 
of two single kinks: either K BB and K AB or K AA and K AB . 

4.3 Kink Variety in Regime II 
Singular kinks 

On the axis we again find the solutions (|69p within the appropriate range of parameters. 

1(a). K^ ±B± : When a < \4>(x)\ < r, the kink orbits connect the A with the B points and 
viceversa. The energy of these kinks is: 

E[K AB ] = \W(±a,0) - W(±t,0)\ = ^ (a + r) 3 (3tt -a 2 - r 2 ) . 

15 

From the arguments of Section §3 one can identify these solitary waves as composite kinks formed 
by two basic kinks. One subtle point is the following: the kink orbits pass through one of the foci. 
Nevertheless, on the <f)± axis the foci are regular points of the superpotential. Although a pencil of 
kink orbits meet at the foci, the K AB trajectories do not belong to this congruence. We conclude 
that these composite kinks are stable. 

1(b). Kf* *: If a < |0i(x)| (|69p . provides another type of solitary waves. The kink orbits 
are also located on the <pi axis but link the B + and B^ points and viceversa, see Figure 11. Their 
energy is: 

E[K? ±B *] = \W(a,0)-W(-a,0)\ = ^a 3 (5r 2 -a 2 ) 
These kink orbits do not cross the foci and are thus stable. 

On the ellipse ([6?]) (left) two types of kinks live, ruled by the formula ([TO]) . 



2(a). K AC : If 1 < |</>i(x)| < r and <f>2{x) = ±G\Jt 2 — 4> 2 (x) the kink orbits connect the A + 
and the C± points (and also the A_ and the C±), see Figure 11. The energy of these kinks is: 

e[k£ c ] = \w(t,o)- w(i,±n)| = |w(-T,o)-w(-i,±n)| 

/2 

= — (T 3 (r-l) 2 (2(j 2 + 4cjO-10T 2 + 6O 2 -5T 3 + 3fi 2 T) 
15 

The kink orbits do not pass through the foci and these kinks are stable. 

2(b). Ki c+ : If 1 < |0i (x) |, formula ()70|) provides the other type of kinks with orbits in the 
quadric curve (|67j) (left) . These kink orbits connect the points C± with CJ, see Figure 11. The 

kink energy is: E[K^' C+ ] = \W(-l,±ft) - W(l,±Cl)\ = ^cr 3 (5r 2 - a 2 ) and these kinks are 
also stable. In this regime the quadric ([67]) (right) is a hyperbola: 

# - t#t = - 2 • ( 74 ) 

3. K± B± : When 1 > |0i(a;)| > o and <p2(%) = ±|f|-\/0 2 (a;) — a 2 , formula ([71]) gives these 
kinks. The kink orbits connect with B_ and C± with B + , see Figure 11. The energy of 
these stable kinks is: E[K? B ] = \ W(-l, ±fi) - W(-a, 0)| = |W(l,±n) - W(a,0)\ = E[K AB ] - 
E[K AC ']. 
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Generic kinks 



In this regime of the model, the kink orbit and the kink profile general equations (|38p and ([39 
provided by the Hamilton-Jacobi method become respectively: 



U — 0~\ t 2 (t 2 -ct 2 ) f U + Vt \ ~S T 7 I I ' T — U \ ct 2 (t 2 -<t 2 ) 
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r + u 



Sign(«') 



(76) 



,2^(^+72 ) 



In M 2 , the boundaries of the cells are formed by pieces of the ellipse (|67|) (left), the hyperbola 
()74p . and the </>i axis, see Figure 15. Although there are three cells in this regime, two of them are 
related by the symmetries of the model. The analysis is thus restricted to the C 23 n and C 2 f cells. 
From the curves determined by (|75p and (|76p . we find the following three types of kink families: 

4(a). A 2 T (7i): In the cell C 23 n the trajectories determined by ([75]) and ([76]) connect the 
points with C_ -skipping the saddle points 1?_ and C+- or B- with -skipping the saddle 
points B + and Cj>, depending on the sign of the combinations Sign(u') = Sign(t/) and Sign(u') ^ 

Sign(u'). The orbits of these K 2 T (7i) one-parametric kink families do not cross the foci, see 
Figure 15, and the kinks are stable. 
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Figure 15: K 2 ± T : (7i) kinks: Profiles and energy density (left) Orbits in R 2 (center) in E 2 (right). 
The energy of any member of these families of topological kinks and the kink mass sum rules are: 



S[< ±C ^( 71 )] 



|W(±<r,0)-W(l,=Fn)| ± ] + E[K^ ± ]=E[k" ± ° ±] + E[K^ b ±] 

/2 

— (-4a 5 + 2(M7 2 + 10O 2 r + 150 3 - 5<r 2 ft 3 - 2r 5 - 50 3 r 2 + 3ft 5 ) 
15 



This result perfectly fits the infinite 7 i limit of the family: 

rB±C) 



hm < ±Ct ( 7 i) 

71— >±co 



± + Af^ c± 

Af ±C±+ 



A^ B± 
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4.(b) K 2 +B+ (^{) and K^ 3 ' (71): Inside the cells Cff and C£|_i, the trajectories determined 
by the equations (|75p and (|76p for any sign combination respectively connect the A + (A_ ) and B + 
(B_) nodal points without crossing any focal point, see Figure 16. These heteroclinic orbits thus 
provide a family of stable topological kinks (one family per cell). 



01 (x) J" 




02 (X^-S- 





-a 
-53 

-n 


X 




Figure 16: -K"^ c (7i) kinks: Profiles and energy density (left) Orbits in the Cartesian plane (center) and 
elliptic plane (right). 
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(a + r) 3 (3Q -a 2 - r 2 



The energy of these topological kinks is: 

E[K^ B± ( 71 )] = \W(±t,Q) - W(±a,0] 
The kink mass sum rule 

E[K^ B± ( 71 )] = E[K? C ] + E[K? B ] = E[K£ B ] 
means that the infinite 71 limit of the family is: 

AC , jsCB 



5. K 



c+f. v^-c: 



(71) and K, 



c+F+yF+c- 



(71 



In the C 2 f and C 2 2_! cells equations ([75]) and 



([75]) include trajectories passing from one C + point to a focus and from a focus to a point ( 
or viceversa) for any combination of signs. Kink orbits (homoclinic in E 2 , heteroclinic in IR 2 ) are 
formed by continuously gluing at a focus two of these trajectories, one with Sign(n)' = Sign(u)' 
and the other with Sign(u)' 7^ Sign(i>)'. In this way, two families of topological kinks are obtained 
with orbits starting from one of the points, crossing the F± focus, and ending at the correlative 
point C±, see Figure 17. The energy of these unstable topological kinks is: 



|w(i,±n) - w(±n,o)| + \w(±n,o)- w(-i,±n)| 

Ay/2, 
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(a + r) d (3fi-r 2 -a 2 



The kink mass sum rule is E[K A 
with the 71 — > 00 limit of the family: 



C+F\J FC~ 



( 7l )] = £[iff B ]+£[j^( 7l )]+£[2Tf G '] in agreement 



"AB 



ACi 



lim KZ +F V FG - M = K? B + Kf B M + Kt 

71— >±oo 



Ac 



In sum, the kink variety in this regime is formed by seven kinds of solitary waves. (1) Four singles 
kinks: Kf 3 , K? b , K BB and K?° '. (2) Three kink families: i^ c (7i), ^ B (7i) and K% c (ji). 
Members of the K BC (^1) family are composite kinks formed by a pair of either K? C -K% B or K BC - 
K BB kinks. The basic kink inter-center distance in each pair is finite and is determined by the 
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etc 



Figure 17: K 4 + + (7i) kinks: Profiles and energy density (left) Orbits inR (center) in E (right). 



71 parameter. Basic kinks appear liberated from each other (infinite inter-center distance) when 
I71I = 00. The kinks of the K^i^i) family are similar non- linear combinations of two single kinks: 
the Kf°-Kf B kink pair arise at the 71 = 00 boundary of the K^iji) moduli space. There is, 
however, a subtle point: at the 71 = —00 boundary lives the K AB kink, which is thus a composite 

C F \/ FC 

but singular kink. Finally, the structure of the K A (71) kinks is even more complex. These 

kinks can be thought of as non-linear combinations of the K AC , K BA and K BC basic kinks when 
the distances between the three centers are finite. 



4.4 Kink Variety in Regime III 
Singular kinks 

In this regime the ansatz <f>i = plugged into (|69p gives rise to two types of kink orbits. 

1(a). K A±B± : These orbits connect the points A + with B + or with 5_, see Figure 12. The 
energy of the corresponding kinks is: E[K^ B ] = \W(±l, 0) - W(±n, 0)| = ^(1 - 0) 3 (1 + 0(3 + 
«))• 

1(b). K 1 ~: There is also a kink orbit on the </>i-axis connecting B + with B-. The kink 
energy is: E[K B+B ~] = \W(to,G) - W{-fl,0)\ = ^0 3 (5 - O 2 ). 

The third type of singular kink in regime III lives on the ellipse (|67|) (left). 



2. K A+A : These elliptic kink orbits are given by ()70l and) connect the A + and A- points, see 



Figure 12. The kink energy is: E[K 1 



|i^(i,o) - w(-i,o)| = ^o 3 (5 - n 2 ) = e\k[ 
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Generic kinks 

The Hamilton-Jacobi procedure provides the kink orbits 
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and the kink profiles 
< 
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of Regime III in the Ci^j cell. 

3. K^~ B+ (71) and iT.f +B ~ (71): Equations ([77|) and (f75|) describe two families of topological kinks 
joining the A± with the -B^ points, see Figure 16.. 




Figure 18: K A B+ (7i) kinks: Profiles and energy density (left) Orbits inM 2 (center) in E 2 (right). 
The kink energy is: 



E[K^( 7l )] = \W(±1,0) - w( T n, 
such that the kink mass sum rules hold: 



2^2 

~l1T 



(1 - 5ft 2 + 15ft 3 - 3ft 5 



E[Kz ±B *( 7l )} = E[K^' A+ ] + E[K A+B+ ] = E[K B ~ B+ \ + E[K A - B -} 
in agreement with the |7i| =00 limit of the family 
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4.5 Kink mass sum rules 

A summary of the different kink mass sum rules is as follows: 

- Kink mass sum rules in regime I: 

E[K AA ] = E[K BB ) 
E[K AB { 11 )]=E[K AB ]+E[K BB ] 
E[K BB ( 7l )] = E[K BFyFF ^ FB ] + E[K BB ] 

- Kink mass sum rules in regime II: 

E[K BB ] = E[KC +C+ ] 

E[K B ^( ll )}=E[K B C + ] + E[K C 1 * 
E[K AB ] = E[K AC+ ] + E[KC +B ] 
E[K AB ( 11 )]=E[K AB ] 

E[K fF\lFC+ {ii)]=2E[K AB {li)] 



CI 
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- Kink mass sum rules in regime III: 

E[K? A ] = E[K BB ] 

EiKfBfri)] = E[K AB ] + E[K BB ] 

5 Summary and Outlook 

The developments discussed in this paper unveil the structure of a broad class of two-component 
scalar field theory models, which generalize the well-studied MSTB model. It is well known that 
the search for static solutions in (1+1)D scalar field theory is tantamount to the solving of an 
analogous mechanical system with flipped potential with respect to the field theoretical model. 
The generalized MSTB models are characterized by: (1) Having a Type I Liouville analogous 
mechanical system. These mechanical systems are Hamilton-Jacobi separable by using elliptic 
coordinates. (2) Having a discrete set of unstable critical points of the mechanical potential, 
homogeneous solutions or "vacua" in the field theory. This feature supports the existence of very 
rich varieties of heteroclinic and homoclinic separatrix trajectories in the mechanical problem, and 
topological and non-topological solitary waves or kinks in the field theory setting. 
We have shown the following structural properties: 

• The "vacua" of the field theoretical system and the foci of the elliptic and hyperbolic coordi- 
nating curves are placed at the junctions of a rectangular mesh displayed on the elliptic strip, 
which is divided into different cells where the dynamics runs separately. 

• The distribution of the vacua on this reticulum depends on the details of the potential energy 
density of the model. 

• There exist only three types of cells in the mesh. In Type I cells, the vacua of the model sit 
on the four vertices. Only three vertices are vacua in Type II cells, whereas one of the two 
foci sits at the fourth vertex. Finally, in Type III cells two vacua and two foci occupy the 
four vertices. 

• The cell of the reticulum is the cornerstone that determines the structure of the kink variety. 
There exist two families of kink solutions confined in each cell and single kinks living on 
the edges of the cell. The characteristics of these kinks, in particular their stability, do not 
depend on the details of the potential density but only on the type of the cell in which they 
are confined. 

• The different kink families are described thoroughly. Because the separability of the mechan- 
ical model, kink mass sum rules arise between at least two types of kinks. 

• For this reason, each member of a kink family (confined inside the cell) is indeed a composite 
kink that can be understood as a non-linear combination of several single kinks (settled at 
the edges of the cell). 

• In sum, the potential energy density of the field theoretical model determines the reticulum 
bounding the cells and the cells determine the behavior and structure of the kink families. 
Therefore, a qualitative description of the kink variety of generalized MSTB models is do-able 
by merely observing the expression of the potential energy density in elliptic coordinates. 

We finish by noticing that the scheme developed in this paper can be extended to other two- 
component scalar field theoretical models. In particular, models with Type III and II Liouville 
analogous mechanical systems, respectively Hamilton-Jacobi separable using parabolic and polar 
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coordinates, have very similar structures. We believe, however, that the Type I models are primor- 
dial because parabolic and polar coordinates are special limits of the system of elliptic coordinates. 
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